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■ ABSTRACT 

^ I Due to their large dynamical mass-to-light ratios, dwarf spheroidal galaxies (dSphs) are 

^SJ ' promising targets for the indirect detection of dark matter (DM) in 7-rays. We examine their 

detectability by present and future 7-ray observatories. The key innovative features of our 
■"nJ" , analysis are; (i) We take into account the angular size of the dSphs; while nearby objects 

■ have higher 7 ray flux, their larger angular extent can make them less attractive targets for 
T^lj- I background-dominated instruments, (ii) We derive DM profiles and the astrophysical J-factor 

. (which parameterises the expected 7-ray flux, independently of the choice of DM particle 

model) for the classical dSphs directly from photometric and kinematic data. We assume very 
little about the DM profile, modelling this as a smooth split-power law distribution, with 
and without sub-clumps, (iii) We use a Markov Chain Monte Carlo (MCMC) technique to 
marginalise over unknown parameters and determine the sensitivity of our derived J-factors to 
both model and measurement uncertainties, (iv) We use simulated DM profiles to demonstrate 
Jh ' that our J-factor determinations recover the correct solution within our quoted uncertainties. 

. . .' Our key findings are: (i) Sub-clumps in the dSphs do not usefully boost the signal; (ii) The 

sensitivity of atmospheric Cherenkov telescopes to dSphs within ^ 20 kpc with cored halos 
can be up to ^ 50 times worse than when estimated assuming them to be point-like. Even for 
the satellite-borne Fermi-LAT the sensitivity is significantly degraded on the relevant angular 
scales for long exposures, hence it is vital to consider the angular extent of the dSphs when 
selecting targets; (iii) No DM profile has been ruled out by current data, but using a prior on 
the inner dark matter cusp slope ^ 7prior ^ 1 provides J-factor estimates accurate to a 
factor of a few if an appropriate angular scale is chosen; (iv) The J-factor is best constrained 
at a critical integration angle etc = 2rhaif /d (where rhaif is the half light radius and d is the 
distance to the dwarf) and we estimate the corresponding sensitivity of 7-ray observatories; 
(v) The 'classical' dSphs can be grouped into three categories: well-constrained and promising 
(Ursa Minor, Sculptor, and Draco), well-constrained but less promising (Carina, Fornax, and 
Leo I), and poorly constrained (Sextans and Leo II); (vi) Observations of classical dSphs with 
Fermi-LAT integrated over the mission lifetime are more promising than observations with 
the planned Cherenkov Telescope Array for DM particle mass < 700 GeV. However, even 
Fermi-LAT will not have sufficient integrated signal from the classical dwarfs to detect DM 
in the 'vanilla' Minimal Supersymmetric Standard Model. Both the Galactic centre and the 
'ultra-faint' dwarfs are likely to be better targets and will be considered in future work. 

Key words: astroparticle physics — (cosmology:) dark matter — Galaxy: kinematics and 
dynamics — 7-rays: general — methods: miscellaneous 
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1 INTRODUCTION 

The detection of 7-rays from dark matter (DM) annihil ation is one 
of the most promisin g channels for indirect detection jOunn et al] 
1 19781 : ISteckeil f 19781). Since the signal goes as the DM density 
squared, the Galactic cent re seems to be the obv ious location to 
search for such a signal dSilk & BloemenI [l987t) . However, it is 
plagued by a confusin g background of astrophysical sources (e.g. 
lAharonian et al .|2004 . For this reason, the dwarf spheroidal galax- 
ies (dSphs) orbiting the Milky Way have been flagged as favoured 
targets given their pot entially high DM densities an d small astro- 
physical backgrounds ( lLakeiri990l : lEvans et al.l2004h . 

Despite the growing amount of kinematic data from the clas- 
sical dSphs, the inner parts of their DM profiles remain poorly 
constrained an d can generally accommodate both cored or cuspy 
soluti ons (e.g. iKoch et alj|2007l : IStrigari et aT]|2007l : IWalker et all 
I2OO9I) . There are two dSphs — Fornax and Ursa Minor — that 
show indirect hints of a cored distribution dKlevna et all l2003l : 
iGoerdt et al.ll2006l) : however, in both cases the presence of a core 
is inferred based on a timing argument that assumes we are not 
catching the dSph at a special moment. Theoretical expectations 
remain similarly uncertain. Cusps are favoured by cosmological 
models th at model the DM alone, assumin g it is cold and collision- 
less (e.g. iNavarro. Frenk & White! IT996I) . However, the complex 
dynamical interplay between stars, gas and DM during galaxy for- 
mation could erase such cusps leading to cored distributions (e.g. 
Navarro et all 1 19 96: Read & Gilmore 2005: Mashchenko et all 
20081 : iGoerdteTaklboiOl : iGovemato et al.ll2010l : ICole et al.ll2oTTir 



Cores could also be an indication of other possibilitie s such 
as self-interacting dark matter (e.g. iHogan & DalcantonI I2OO0I : 
iMoore et al.ll2000l) . 

Knowledge of the inner slope of the DM profile is of crit- 
ical importance as most of the annihilation flux comes from 
that region. Lacking this information, several studies have fo- 
cused on the detectability of these dSphs by cuixent 7-ray ob- 
servatories such as the satellite-borne Fermi-LAT and atmo- 
spheric Cherenkov telescopes (ACTs) such as H.E.S.S., MAGIC 
and VERITAS, using a small sample of cusped and cored pro- 
files (generally one of each). Most studies rely on standard core 
and cus p profiles fitted to the kin e matic data of the dSph of 
interest 1 Bergst rom & H oopei 20061: Sanchez-Conde et"an |2007| : 
iBringmann 673^.2009, : .Pieri et al.ll2009l : iPieri et al.ll2009l) . Other 
authors use a 'co smological pri or' from large scale cosmological 
simulations (e.g. |Kuhlen' 2010). Both approaches may b e com - 
bined, such as in [strigari et al. (2007) and Martin ez et al. I (l200^ 
who rely partially on the results of structure formation simulations 
to constrain the inner slope and then perform a fit to the data to 
derive the other parameters. However such cosmological priors re- 
main sufficiently uncertain that their use is inappropriate for guid- 
ing o bservational strat egies. There have been only a few studies 
(e.g., lEssig et"Zll2009l) which have not assumed strong priors for 
the DM profiles. 

In this work, we revisit the question of the detectability of dark 
matter annihilation in the classical Milky Way dSphs, motivated by 
ambitious plans for next-generation ACTs such as the Cherenkov 
Telescope Array (CTA). We rely solely on published kinematic data 
to derive the properties of the dSphs, making minimal assumptions 
about the underlying DM distribution. Most importantly, we do not 
restrict our survey of DM profiles to those suggested by cosmolog- 
ical simulations. We also consider the effect of the spatial extent of 
the dSphs, which becomes important for nearby systems observed 
by background-limited instruments such as ACTs. 



This paper extends the earlier study of IWalker et al. I teoiil) 
which showed that there is a critical integration angle (twice the 
half-light radius divided by the dSph distance) where we can ob- 
tain a robust estimate of the J-factor (that parameterises the ex- 
pected 7-ray flux from a dSph independently of the choice of dark 
matter particle model; see Section |2j, regardless of the value of 
the central DM cusp slope 7. Here, we focus on the full radial de- 
pendence of the J-factor. We consider the effect of DM sub-lumps 
within the dSphs, discuss which dSphs are the best candidates for 
an observing programme, and examine the competitiveness of next- 
generation ACTs as dark matter probes. 

This paper is organised as follows. In Section |2] we present a 
study of the annihilation 7-ray flux, focusing on which parameters 
critically affect the expected signal. In Section |3] we discuss the 
sensitivity of present/future 7-ray observatories. In Section |4] we 
present our method for the dynamical modelling of the observed 
kinematics of stars in dSphs. In Section |5] we derive DM den- 
sity profiles for the classical dSphs using an MCMC analysis, from 
which the detection potential of future 7-ray observatories can be 
assessed. We present our conclusions in Section |60 

This paper includes detailed analyses from both high-energy 
astrophysics and stellar dynamical modelling. To assist readers 
from these different fields in navigating the key sections, we sug- 
gest that those who are primarily interested in the high-energy cal- 
culations may wish to focus their attention on Sections |2]|3] and |5] 
before moving to the conclusions. Readers from the dynamics com- 
munity may instead prefer to read Sections|2]|4]and|5] Finally, those 
who are willing to trust the underlying modelling should proceed 
to Section [5] where our main results regarding the detectability of 
dSphs are presented in Figs. [12] [15] ll6land[T7l 



2 THE DARK MATTER ANNIHILATION SIGNAL : KEY 
PARAMETERS 

2.1 The 7-ray flux 

The 7-ray flux <1?^ (photons cm^^ s^^ GeV^'^) from DM annihi- 
lations in a dSph, as seen within a solid angle Afi, is given by (see 
Appendix IaI for definitions and conventions used in the literature): 



^(iJ^,An)=$PP(i5,)x J(An), 



(1) 



The first factor encodes the (unknown) particle physics of DM an- 
nihilations which we wish to measure. The second factor encodes 
the astrophysics viz. the l.o.s. integral of the DM density-squared 
over solid angle ASl in the dSph — this is called the 'J-factor'. We 
now discuss each factor in turn. 



^ Technical details are defeiTed to Appendices. In Appendix IaI we com- 
ment on the various notations used in similai' studies and provide conversion 
factors to help compare results. In Appendix|B] we provide a toy model for 
quick estimates of the J-factor. In Appendix|D] we calculate in a more sys- 
tematic fashion the range of the possible 'boost factor' (due to DM clumps 
within the dSphs) for generic dSphs. In Appendix [E] we show that con- 
volving the signal by the PSF of the instrument is equivalent to a cruder 
quadrature sum approximation. In Appendix |F] we discuss some technical 
issues related to confidence level detennination from the MCMC analysis. 
In Appendix|G] the reconstraction method is validated on simulated dSphs. 
In Appendix [H] we discuss the impact of the choice of the binning of the 
stars and of the shape of the light profile on the J-factor determination. 
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2. 1. 1 The particle physics factor 

The particle physics factor is given by: 



2mJ dE^ 



(2) 



where is the mass of the DM particle, Uann is its self- 
annihilation cross-section and (uannw) the average over its velocity 
distribution, and dA'^^ /d-E^ is the differential photon yield per an- 
nihilation. A benchm ark value is (cTann^') ~ 3 x 10~^^ cm"^ s^^ 
jjungman et al.|[l996h . which would result in a present-day DM 
abundance satisfying cosmological constraints. 

Unlike the annihilation cross section and particle mass, the 
differential annihilation spectrum (dA'^/d_E-y(_E^)) requires us to 
adopt a specific DM particle model. We focus on a well-motivated 
class of models that are within reach of up-coming direct and in- 
direct experiments: the Minimal Supersymmetric Standard Model 
(MSSM). In this framework, the neutralino is typically the light- 
est stable particle and therefore one o f the most favoured DM can- 
didates (see e.g. ISertone et al.ll200^ . A 7-ray continuum is pro- 
duced from the decay of hadrons (e.g. tv" — >■ 77) resulting from 
the DM annihilation. Neutralino annihilations can also directly pro- 
duce mono-energetic 7-ray lines through loop processes, with the 
forma tion of either a pair of 7-rays (xx ^ 77- Bergstr om & UUiol 
1997 ). or a boson and a 7-ray (xx ~^ 'YZ°: IUllio & BergstromI 
19981) . We do not take into account such line production processes 
since they are usually sub-dominant and very model dependent 
teringmann et al.l|2008l) . The differential photon spectrum we use 
is restricted to the continuum contribution and is written as: 



{E,)^J2^, 



dm, 

dE. 



(3) 



where the different annihilation final states i are characterised by a 
branching ratio 6i. 

Using the parameters in iFomengo et al.l l l2004l) . we plot the 
continuum spectra calculated for a 1 TeV mass neutralino in Fig. 
[T] Apart from the r+r^ channel (dash-dotted line), all the an- 
nihilation channels in the continuum result in very similar spec- 
tra of 7-rays (dashed lines). For charged annihilation products, 
internal bremsstrahlung (IB) has recently been investigated and 
found to enhance the s pectrum close to the kinematic cut-off (e.g., 
iBringmann et aUFlOOSh . As an illustration, the long-dashed line in 
Fig. [T] corresponds to the benchmark configuratio n for a wino- 
like neutralino taken from IBringmann et al.l fcOOSh . However, the 
shape and amplitude of t his spectrum are str ongly model de p enden t 
teringmann et al.l |2009|) and, as argued in ICannoni et al.l ( |2010|) . 
this contribution is relevant only for models (and at energies) where 
the line contribution is dominant over the secondary photons. 

We wish to be as model-independent as possible, and so do not 
consider internal bremsstrahlung. In the remainder of this paper, all 
our results will be based on an averas e spectrum taken from the 
parametrisation teergstrom et al.|[l998l solid line in Fig.[T): 



diV-, 
dE^ 



^_dNj_ 

m-y dx 



1 0.73 e" 



(4) 



with X = E-y/m-^. Finally, in order to be conservative in deriving 
detection limits, we also do not consider the possible 'Sommerfeld 
enhancement' of the DM annihilation cross-section l lHisano et al.l 



' ' ' ' ' ' ' I 
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Figure 1. Differential spectra (multiplied by x^) of 7-rays from the frag- 
mentation of neutrino annihilation products (here for a DM particle mass 



iFomengo et al.1 b004) and an average parametrisation iBergstrom et alj 
119981) is marked by the black solid line; this is what we adopt through- 
out this paper. The black dashed Hne is the benchmark model BM4 
iBringmann et alj20 Q8) which includes internal bremsstrahlung and serves 
to illustrate that veiy different spectra are possible. However, the example 
shown here is dominated by line emission and therefore highly model de- 
pendent; for this reason, we do not consider such effects in this paper 



l2004l2005l) PI This depends inversely on the DM particle velocity, 
and thus requires precise modelling of the velocity distribution of 
the DM within the dSph; we will investigate this in a separate study. 

2.1.2 The J-f actor 

The second term in Eq. ^ is the astrophysical J-factor which de- 
pends on the spatial distribution of DM as well as on the beam size. 
It corresponds to the l.o.s. integration of the DM density squared 
over solid angle Afi in the dSph: 



J- 



Puuil-, ^) dldO,. 



(5) 



The solid angle is simply related to the integration angle a^^t by 

AQ = 27r ■ (1 - COs(Qint)) . 

The J-factor is useful because it allows us to rank the dSphs by 
their expected 7-ray flux, independently of any assumed DM par- 
ticle physics model. Moreover, the knowledge of the relative J- 
factors would also help us to evaluate the validity of any poten- 
tial detection of a given dSph, because for a given particle physics 
model we could then scale the signal to what we should expect to 
see in the other dSphs. 

All calculations of J presented in this paper were performed 
using the publicly available CLUMPY package (Charbonnier, Com- 
bet, Maurin, in preparation) which includes models for a smooth 
DM density profile for the dSph, clumpy dark matter sub-structures 
inside the dSph, and a smooth and clumpy Galactic DM distribu- 
tion.0 

^ This effect depends on the mass and the velocity of the particle; the re- 
sulting boost of the sig nal and the impac t on detectability of the dSphs has 
been discussed, e.g., in IPieri et aljf2009l) . 

In Appendix|B] we provide approximate formulae for quick estimates of 
the J-factor and cross-checks with the numerical results. 
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2.1.3 DM profiles 

For the DM halo we use a generalised (a, /3, 7) He mquist profile 
given by l lHemquisJll990l : lOehnenll 19931 : IZhadI 19961) : 



p{r) = pA — 



l+(^ 



(6) 



where the parameter a controls the sharpness of the transition 
from inner slope, limr^o dln(p)/dlii(r-) = —7, to outer slope 
linir^oo dln(p)/dln(r) = — /3, and r-g is a characteristic scale. 
In principle we could add an additional parameter in order to in- 
troduce an exponential cut-off in the profile of Eq (|6) to mimic 
the effects of t idal truncation, as prop osed in e.g., the Aqu arius 
dSpringel et"al]|2008..) or Via Lactea II jPiemand etal]|2008l) sim- 
ulations. However, the freedom to vary parameters Vs, a and P 
in Eq ([6f already allows for density profiles that fall arbitrarily 
steeply at large radius. Moreover, given that our MCMC analysis 
later shows that the outer slope /3 is unconstrained by the available 
data and that the J-factor does not correlate with /3, we choose not 
to add further shape parameters. 

For profiles such as 7 1.5, the quantity J from the inner 
regions diverges. This can be avoided by introducing a saturation 
scale Tsat , that corresponds physically to the typical scale where the 
annihilation rate [(o^^')p(rsat)/w^y] balances the gravita tional 
infall rate of DM particles (Gp)"^/^ jBerezinskv et aljl992h . Tak- 
ing p to be about 200 times the critical density gives 

'10"^®cm^ s 



Psat 



3 X 10" 



100 GeV 



{av) 



Mq kpc 



The associated saturation radius is given by 

1/7 



^sat 



Ps 
Psat 



(7) 



(8) 



This limit is used for all of our calculations. 



2.2 Motivation for a generic approach and reference models 

In many studies, the 7-ray flux (from DM annihilations) 
is calculated using the point-so u rce a pproximation (e.g., 
iBergstrom & Hooperl I2OO6I: iKuhlerJ I2OI0I) . This is valid so 
long as the inner profile is steep, in which case the total luminosity 
of the dSph is dominated by a very small central region. However, 
if the profile is shallow and/or the dSph is nearby, the effective 
size of the dSph on the sky is larger than the point spread function 
(PSF) of the detector, and the point-source approximation breaks 
down. For upcoming instruments and particularly shallow DM 
profiles, the effective size of the dSph may even be comparable 
to the field of view of the instrument. This difference in the 
radial extent of the signal does matter in terms of detection (see 
Section[3](. Hence we do not assume that the dSph is a point-source 
but rather derive sky-maps for the expected 7-ray flux. 



2.2.1 Illustration: a cored vs ciisped profile 

Fig. |2] shows J as a function of the integration angle aint for a 
dSph at 20 kpc (looking towards its centre). The black solid line is 
for a cored profile (7 — 0) and the green dashed line is for a cuspy 
profile(7 — 1.5); both are normalised to unity at Qfint ~ 5°.Forthe 
cuspy profile, ~ 100% of the signal is in the first bin while for the 
cored profile, J builds up slowly with Oint, and 80% of the signal 
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Figure 2. Finite size effects: J as a function of the integration angle ajnt 
for a dSph at 20 kpc (pointing towards the centre of the dSph). The black 
solid line is for a cored profile (7 = 0) and the green dashed line is for a 
cuspy profile (7 = 1.5): both are nomialised to unity at a-^t = 5°. 



also indicated by the symbols which show the contribution of DM 
shells in two angular bins — whereas the (green) hollow squares 
have a spiky distribution in the first bin (7 = 1.5), the (black) filled 
circles (7 = 0) show a very broad distribution for J. 

The integration angle required to have a sizeable fraction of 
the signal depends on several parameters: the distance d of the 
dSph, the inner profile slope 7, and the scale radius Vs. Small in- 
tegration angles are desirable since this minimises contaminating 
background 7-ray photons and maximises the signal to noise. Thus 
the true detectability of a dSph will depend on its spatial extent on 
the sky, and thus also on d, 7 and rs. 

2.2.2 Generic dSph profiles 

As will be seen in Section[5] the errors on the density profiles of the 
Milky Way dSphs are large, making it difficult to disentangle the 
interplay between the key parameters for detectability. Hence we 
select some 'generic profiles' to illustrate the key dependencies. 

The most constrained quantity is the mass within the half-light 
radius rhaif (typically a few tenths o f a kpc), as this is where most of 
the ki nematic data come from (e.g.. I Walker et al ]|2009l:IWolfetal.l 
I2OIOI) . For the classical Milky Way dSphs, the typical mass within 
fhaif ~ 300 pc is found to be M300 ~ IO'^Mq jStrigari et al.l2008l 
— see also the bottom panel of Fig. I13t . If the DM scale radius 
is significantly larger than this (rs ^ rhaif) and the inner slope 
7 > 0.5, we can approximate the enclosed mass by: 



M3, 



47r psrg / 300pc 
3 - 7 V ''s 



3-7 



lO^M, 



■ 



(9) 



(w.r.t. the value for aint = 5°) is obtained for a' 



80% 



3°. This is 



The parameter ps is thus determined completely by the above con- 
dition, if we choose the scale radius rs and cusp slope 7. 

Table [T] shows, for several values of rs and 7, the value re- 
quired for Ps to obtain the assumed A/300 mass. We fix a = 1, /3 = 
3 but our results are not sensitive to these choices. QThe values of rs 
are chosen to encompass the range of rs found in the MCMC anal- 
ysis (see Section[5)- To further convince ourselves that the generic 



^ For a different mass for the dSph, the results for J below have to be 
reseated by a factor (Alggo /W'^Mq)'^ since the density is proportionat to 
M300 , while J goes as the density squared. 



Dark matter annihilation in dwarf spheroidal galaxies and j-ray observatories: I. Classical dSphs 5 



Table 1. The required normalisation ps to have Afsoo = IO^Mq for a 
sample of (1, 3, 7) profiles with varying scale radius Tg. 



Pa {WMq kpc-3) 
7 \ rs [Icpc] 0.10 0.50 1.0 
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88 


3.2 
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profiles we present here are a possible description of real dSphs, we 
checked (not shown) using typical stellar profiles and properties of 
these objects (i.e., halflight radius of a few 100 pc), that a flat ~ 10 
km s^^ velocity dispersion profile within the eiTor bars is recov- 
ered. We also study below the effect of moving these dSphs from 
a distance of lOkpc to 300 kpc, corresponding to the typical range 
covered by these objects. 

2.2.3 Sub-structures within the dSph 

Structure formation simulations in the currently favoured ACDM 
(cold DM plus a cosmological constant) cosmology find that DM 
halos are self-similar, containing a we alth of smaller 'sub-str ucture' 
halos down to Earth-mass halos (e.g. iDiemand et al. I I200I . How- 
ever, as emphasised in the introduction, such simulations typically 
neglect the influence of the baryonic matter during galaxy forma- 
tion. It is not clear what effect these have on the DM sub-structure 
distribution. For this reason, we adopt a more generic approach. We 
assess the importance of clumps using the following recipejf] 

(i) we take a fraction / = 20% of DM mass in the form of 
clumps; 

(ii) the spatial distribution of clumps follo ws the smooth one; 

(iii) the clump profiles are calculated a /fliBullock etaljj200lh 
(here after BOl), i.e. an 'NFW' profile dNavarro, Frenk & White! 
[199^1 with concentration related to the mass of the clumps. 

(iv) the clump mass distribution is oc M"" (a = —1.9), within 
a mass range M^in - Af^ax = [10"*^ - 10*^] Mq. 

Although these parameters are very uncertain, they allow us 
to investigate the impact of substructures on the J-factor. They are 
varied within reasonable bounds in S ection |2 . 3 . 2 K and AppendixiDt 
to determine whether the sub-clump contribution can boost the 
signal. Note that a 20% clump mass fraction is about twice as 
large as the fraction o btained from numerical simulations (see, e.g., 
ISpringeletal]|2008l) . This generous fraction does not affect our 
conclusions, as discussed below. 

2.3 Jsm and Jsubci for the generic models 

As an illustration, we show in Fig.|3]one realisation of the 2D dis- 
tribution of J from a generic core profile (7 = 0) with 7-3 = 1 
kpc (sub-clump parameters are as described in Section |2.2.3| l. The 

^ More details about the clump distributions can be found in Appendix lB2l 
See also, e.g.. Section 2 in Lavalle et al. (2008) and references therein, as 
we use the same definitions as those given in that paper. 



dSph is at rf — 100 kpc. We note that our consideration of a 
7 = smooth component with NFW sub-clumps is plausible if, 
e.g., baryon-dynamical processes erase cusps in the smooth halo 
but cannot do so in the sub-subhalos. The total J is the sum of the 
smooth and sub-clump distributions. The centre is dominated by the 
smooth component, whereas some graininess appears in the out- 
sldrts of the dSph. In this particular configuration, the 'extended' 
signal from the core profile, when integrated over a very small 
solid angle, could be sub-dominant compared with the signal of 
NFW sub-clumps that it hosts. The discussion of cross-constraints 
between detectability of sub-halos of the Galaxy vs. sub-clumps in 
the dSph is left for a future study. 

In the remainder of the paper, we will replace for simplicity 
the calculation of Jsubci (ctint) by its mean value, as we are pri- 
marily interested in 'unresolved' observations. Hence clumps are 
not drawn from their distribution function, but rather (Jsubci) is 
calculated from the integration of the spatial and luminosity (as a 
function of the mass) distributions (see Appendix IB2t . 

2.3.1 Radial dependence J{9) 

The radial dependence of J is shown in Fig. |4] for four values of 
7 (for an integration angle aint ~ 0.01°). The dashed lines show 
the result for the smooth distribution; the dotted lines show the sub- 
clump contribution; and the solid lines are the sum of the two. The 
peak of the signal is towards the dSph centre. As long as the dis- 
tribution of clumps is assumed to follow the smooth one, regard- 
less of the value of 7, the quantity (1 — /)^ Jsm(O) always domi- 
nates (at least by a factor of a few) over ( Jsubci(O)). (Recall that in 
our generic models, all dSphs have the same A/goo.) The scatter in 
Jtot (0) is about 4 orders of magnitude for 7 G [0.0 — 1.5], but only 
a factor of 20 for 7 £ [0.0 — 1.0]. Beyond a few tenths of degrees, 
(Jsubci) dominates. The crossing point depends on a combination 
of the clump mass fraction /, 7, rs, d, aint. The dependence of J 
on the two latter parameters are discussed in Appendix |C] The ra- 
dial dependence is as expected: the smooth contribution decreases 
faster than that of the sub-clump one, because the signal is propor- 
tional to the squared spatial distribution in the first case, but directly 
proportional to the spatial distribution in the second case. Halving / 
to match the fraction from N-body simulations would have a 25% 
effect on (1 — Jam, but decrease Jsubci by a factor 4, so that 
the cross-over between the two components would occur at a larger 
angle in Fig.|4l 

2.3.2 Boost factor 

Whether or not the signal is boosted by the sub-clump pop ulation 
is still debated in the literature (iStrigar i et al. 2007; Kuhl en et al.l 
l2008l ; |Pieri e7aLll2008l ; Fieri et alfcoogi) . As underiined in the pre- 
vious sections, the sub-clump contribution towards the dSph centre 
never dominates over the smooth one if the spatial profile of the 
sub-clumps follows that of the smooth distribution, and if the inte- 
gration angle remains below some critical angle discussed below. 

Let us first define properly the parameters with respect to 
which this boost is calculated, as there is sometimes some con- 
fusion about this. Here, we define it with respect to the integration 
angle aint (the pointing direction is still towards the dSph centre): 

R/ ^ _ (1 - /)^Jsm(Qint) + Jsubci (Qint) 

B[aint) = J — 7 T . (10) 

Jsm (aint J 

In most studies, the boost has been calculated by integrating out to 
the clump boundary (i.e., af^t ~ Rvii/d). But the boost depends 
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Figure 3. 2D view {x and y axis are in degrees) of J for the generic dSph 
witli 7 = and rs = 1 kpc at d = 100 kpc (M300 = W Mq). 
The sub-clumps are drawn from the reference model described in Sec- 
tion 12.2.31 i.e. f=20%, sub-clump distribution follows smooth, and sub- 
clump inner profiles have NFW with BOl concentration. From top to bot- 
tom panel: aint = 0.1°, 0.05°, and 0.01°. For the sake of comparison, the 
same colour scale is taken for the three integration angles (J is in units of 
kpc~5). 



10'^ 
10^ 



10° n 




from dSph centi-e [deg] 



Figure 4. J as a function of the angle d away from the dSph centre for a 
dSph at 100 kpc with = 0.5 kpc (ps is given in Table[TJ. The integration 
angle is ajnt = 0.01°. For the four inner slope values 7, the various con- 
tributions to J are shown as solid (total), dashed (smooth), and dotted lines 
(sub-clumps). 
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Figure 5. Boost factor as a function of dint x (d/100 kpc) for profiles 
sub-clumps /b/Zon: smooth (see Section [2.2.3) : the dSph is at d = 100 kpc 
(lines) or d = 10 kpc (symbols). 



crucially on aint (the radial dependence of the smooth and sub- 
clump contributions differ, see Section |2.3.1) . 

We plot in Fig.[5]the boost for different inner slopes 7, where 
a direct consequence of Eq. CIll is the Qint X d rescaling. For 
Ts < 0.1 kpc (regardless of 7), or for 7 > 1.5 (regardless of r^), 
the signal is never boosted. Q For small enough aint, B is smaller 
than unity, and if 7 is steep enough, B ~ (1 — /)^. For large 
values, a plateau is reached as soon as aintd Rvii (taken to be 3 
kpc here). In between, the value of the boost depends on Vs and 7 of 
the smooth component. Going beyond this qualitative description is 
difficult, as the toy model formulae of Appendix IB2I gives results 



^ The difference between the level of boost observed for rs = 0.1 kpc or 
Ts = 1 kpc can be understood if we recall that the total mass of the clump 
is fixed at 300 pc, regardless of the value of 7 or Tg. For rs = 0.1 kpc, pa ~ 
©(IO^Mq kpc-3), whereas for = 1 kpc, ps ~ 0{W Mq kpc'^). 
As Jsm OH Pa whereas J^ub oc Ps , the relative amount of J^ub with respect 
to Jsni is expected to decrease with smaller rg. This is indeed what we 
observe in the figure (solid vs dashed lines). 
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correct to only a factor of ~ 2 (which is inadequate to evaluate the 
boost properly). 

To conclude, the maximum value for sub-clump follows 
smooth is < 2, and this value is reached only when integrating the 
signal out to Rvh/d. The boost could still be increased by varying 
the sub-clump properties (e.g., taking a higher concentration). Con- 
versely, if dynamical friction has caused the sub-clump population 
to become much more centrally concentrated than the smooth com- 
ponent, then the boost is decreased. This is detailed in AppendixiDl 
For the most realistic configurations, there is no significant boost 
when a clump mass fraction f = 20% is used. Naturally this result 
is even more true for the smaller / found in N-body simulations so 
we disregard the boost for the rest of this paper and consider only 
the smooth contribution. 



3 SENSITIVITY OF PRESENT/FUTURE 7-RAY 
OBSERVATORIES 

Major new ground-based 7-ray observatori e s are in the plan- 
ning stage, with CTA dCTA ConsortlumI 1201 Oh and AGIS 
jAGIS Collaboration! I2OI0I) as the main concepts. As the designs 
of these instruments are still evolving, we adopt here generic per- 
foimance curves (described below), close to the stated goals of 
these projects. For the Large Area Telescope (LAT) of the Fermi 7- 
ray satellite, the performance for 1 year observations of point-like, 
high Galactic latitude sources is known (Fermi-LAT Collaboration 
|2010|) . but no information is yet available for longer exposures or 
for extended objects. We therefore adopt a toy likelihood-based 
model for the Fermi sensitivity, tuned to reproduce the 1 year point- 
source curves. We note that whilst this approach results in approx- 
imate performance curves for both the ground- and space-based 
instruments, it captures the key differences (in particular the differ- 
ences in collection area and angular resolution) and illustrates the 
advantages and limitations of the two instrument types, as well as 
the prospects for the discovery of DM annihilation in dSphs within 
the next decade. 

3.1 Detector models 

The sensitivity of a major future 7-ray observatory based on an ar- 
ray of Cherenkov Telescopes (FCA in the following, for 'Future 
Cherenkov Array') is approximated based on the point-source dif- 
ferential sensitivity c urve (for a 5<j detectio n in 50 hours of obser- 
vations) presented bv lBemlohr et al.l(l2008h . Under the assumption 
that th e angular resolutio n of such a detector is a factor 2 better than 
HESS jFunk et alj2008l) and has the same energy-dependence, and 
that the effective collection area for 7-rays grows from lO'^m^ at 
30 GeV to 1 km^ at 1 TeV, the implied cosmic-ray (hadron and 
electron) background rate per square degree can be inferred and 
the sensitivity thus adapted to different observation times, spectral 
shapes and source extensions. Given that the design of instruments 
such as CTA are not yet fixed, we consider that such a simplified 
response, characterised by the following functions is a useful tool 
to explore the capabilities of a generic next-generation instrument: 

LS = -13.1 -0.33X + 0.72Js^:^ (11) 
LA = 6 + 0.46X-0.56X^ (12) 
^68 = 0.038 + exp-(X + 2.9)/0.61, (13) 
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Figure 6. The cone angle encompassing 80% of the annihilation flux at as 
function of the inner slope 7. Several different values of rs and distance d 
are shown for each 7, all scaled by (1 kpc/r's and 100 kpc/d). The best-fit 
curve is also shown, corresponding to Eq. l ll6t . 



where 

X = logio (PhotonEnergy/TeV), (14) 

LS = logio (Differential Sensitivity/erg cm^^ s^^), LA = 
logi()(EffectiveArea/m^), and ipeg is the 68% containment 
radius of the point-spread-function (PSF) in degrees. 

For the Fermi detector a similar simplified approach 
is taken, the numbers used b elow being those provided by 
iFermi-LAT Collaboration! (l2010h . The effective area changes as a 
function of energy and incident angle to the detector, reaching a 
maximum of « 8000 cm^. The effective time-averaged area is then 
eAQ/in and the data-taking efficiency e ~ 0.8 (due to instrument 
dead-time and passages through the South Atlantic Anomaly). The 
point spread function again varies as a function of energy (with a 
much smaller dependence as a function of incidence angle) varying 
from 10 degrees to a few tenths of a degree over the LAT energy 
range. A rate of 1.5 x 10"''' cm~^ s~^ sr"^ (> 100 MeV) and 
a photon index of 2.1 are assumed for the background. The sensi- 
tivity is then estimated using a simplified likelihood method which 
provides results within 20% of the sen sitivity for a one year ob- 
servation of a point-like source given bv lFermi-LAT Collaboration! 

Whilst both detector responses are approximate, the compar- 
ison is still useful. Our work incoiporates several key aspects not 
considered in earlier studies, including the strong energy depen- 
dence of the angular resolution of both ground and space based 
instruments in the relevant energy range of 1 GeV to 1 TeV and 
hence the energy-dependent impact of the angular size of the target 
region. 



3.2 Relative performance for generic hales 

Using the results from Section !2. 2. 2| and the detector performance 
models defined above we can begin to investigate the sensitivity of 
future ACT arrays and the Fermi-LAT detector (over long obser- 
vation times) to DM annihilation in dSphs. The detectability of a 
source depends primarily on its flux, but also on its angular extent. 
The impact of source extension on detectability is dealt with ap- 
proximately (in each energy bin independently) by assuming that 
the opening angle of a cone which incorporates 80% of the signal 
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is given by 



(15) 



where i/)go = 1.25 V'es is assumed for the FCA and interpolated 
from values given for 68% and 95% containment for the LAT 
iFermi-LAT Collaborationl(l2010l) : here ago is the 80% containment 
angle of the halo emission. The validity of this approximation (at 
the level of a few percent) has been tested (see Appendix ^ by 
convolving realis tic halo profi les with a double Gaussian PSF as 
found for HESS ( iHomsllzOOSi) . An 80% integration circle is close 
to optimum for a Gaussian source on a flat background (in the back- 
ground limited regime). Fig. |6] shows the 80% containment radius 
of the annihilation flux of generic halos as a function of the inner 
slope 7. This result can be parametrised as: 



ago = 0.8° (1-0.487-0.1377^) 



d 



Ikpc/ VlOOkpc 



(16) 



It is clear that for a broad range of d, 7 and the characteristic 
angular size of the emission region is larger than the angular reso- 
lution of the instruments under consideration. It is therefore critical 
to assess the performance as a function of the angular size of the 
dSph as well as the mass of the annihilating particle. 

Fig. |7] shows the relative sensitivity of Fermi and an FCA 
within our framework as a function of the mass of the annihilat- 
ing particle, adopting the annihilation spectrum given in Eq. 
with the several panels illustrating different points. From Fig.|7]top 
(the case of a point-like signal for different observation times) it is 
clear that Fermi-LAT has a considerable advantage for lower mass 
DM particles (m^ <C 1 TeV) on the timescale for construction of 
an FCA (i.e. over a 5-10 year mission lifetime) in comparison to 
a deep ACT observation of 200 hours. Furthermore, Fermi-LAT is 
less adversely affected by the angular extent of the target regions 
(see Fig. |7] bottom), due to its modest angular resolution in the 
energy range where it is limited by background, meaning that the 
source extension is well matched to the PSF of the instrument. The 
middle panel of this figure illustrates the impact of different ap- 
proaches to the analysis. In the case that there is a DM candidate 
inferred from the discovery of supersymmetry at the LHC (quite 
possible on the relevant timescale) a search optimised on an as- 
sumed mass and spectral shape can be made (solid curves). How- 
ever, all instruments are less sensitive when a generic search is un- 
dertaken. Simple analyses using all the photon flux above a fixed 
energy threshold (arbitrarily set to reduce background) are effec- 
tive only in a relatively nareow range of particle mass. For example 
keeping only >I00 GeV photons works well for ACTs for 0.3-3 
TeV particles; whereas keeping all photons >1 GeV works moder- 
ately well in the 0.1-0.2 TeV range, but is much less sensitive than 
the higher threshold cut over the rest of the candidate dark matter 
particle mass range. The features of these curves are dictated by 
the expected shape of the annihilation spectrum. From Eq. (|4j the 
peak photon output (adopting the average spectrum for DM anni- 
hilation) occurs at an energy which is an order of magnitude below 
the particle mass - effective detection requires that this peak occurs 
within (or close to) the energy range of the instrument concerned. 

The total annihilation flux from a dSph increases at smaller 
distances as l/d^ for fixed halo mass, making nearby dSphs at- 
tractive for DM detection. However, as Fig.|7]shows, the increased 
angular size of such nearby sources raises the required detection 
flux. Fig. [8] illustrates the reduction in sensitivity for an FCA with 
respect to a point-like source for generic dSph halos as a function 
of distance, for inner slopes, 7, of zero and one and with Ts fixed to 




m, [TeV] 



Figure 7. Approximate sensitivities of Fermi-LAT (blue lines), HESS 
(black lines) and the FCA described above (red lines) to a generic halo 
with J = 10^'^ Mg kpc~^, as a function of the mass of the annihilating 
particle and for the annihilation spectrum of Eq. (|4j. Top: The impact of 
observation time is illustrated: dashed lines give the 1 year and 20 hour sen- 
sitivities for Fenni and FCA/HESS respectively while the solid lines refer 
to 10 year (200 hour) observations. Middle: the impact of analysis meth- 
ods is considered for 5 year (100 hour) observations using Fermi (FCA). 
Sohd hues show likelihood analyses in which the mass and spectrum of the 
annihilating paiticle are known in advance, while dashed and dotted lines 
show simple integral flux measurements above fixed thresholds of 1 GeV 
(dashed) and 100 GeV (dotted). Note that the 1 GeV cut implies accepting 
all events for the FCA (where the trigger threshold is Ri20 GeV). Bottom: 
the impact of the angular extension of target sources, as given by the halo 
profile in Fig.|6]is illustrated. The solid lines reproduce the likehhood case 
from the middle panel for a point-like source, and with values of ago of 
0.1° (dashed) and 1 ° also shown. 
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Figure 8. Relative DM annihilation detection sensitivity for a 100 hourFCA 
observation, as a function of dSph distance for different inner slopes 7 and 
with r's fixed to 1 kpc. The sensitivity for a realistic approach using 9go 
is given relative to the sensitivity to a point-like source with the same flux. 
Larger values correspond to poorer performance (larger values of the min- 
imum detectable flux). The assumed spectral shape is again as given by 
Eq. (4) with m-^ = 300 GeV. This sensitivity ratio depends on the strat- 
egy used to estimate the background level at the dSph position. The dashed 
lines show the impact of using an annulus between 3.5° and 4.0° of the 
dSph centre as a background control region. The solid line assumes that the 
background control region lies completely outside the region of emission 
from the dSph. 



1 kpc, relative to the assumption of the full annihilation signal and 
a point-like source. Even for 7 = 1, the point-like approximation 
leads to an order of magnitude overestimate of the detection sensi- 
tivity for nearby (~ 20 kpc) dSphs. A further complication is how to 
establish the level of background emission arising from the residual 
non-7-ray background. A common method in ground-based 7-ray 
astronomy is to estimate t his background fro m an annulus around 
the target source (see, e.g., lBerge et al.ll2007h . The dashed lines in 
Fig- m show the impact of estimating the background using an an- 
nulus between 3.5° and 4.0° from the target. This approach has a 
modest impact on sensitivity and is ignored in the following dis- 
cussions as it reduces the detectable flux but also 6so and leads to a 
small improvement in some cases. 



4 JEANS/MCMC ANALYSIS OF DSPH KINEMATICS 

4.1 dSph kinematics with the spherical Jeans equation 

Extensive kinematic surveys of the stellar components of dSphs 
have shown that these systems have negligible rotational sup- 
port (with the pos sible exception of the Sculptor dSph, see 
iBattaglia et alj2008l) . If we assume that the dSphs are in virial equi- 
librium, then their internal gravitational potentials balance the ran- 
dom motions of their stars. In order to estimate dSph masses, we 
consider here the behaviour of dSph stellar velocity dispersion as 
a function of distance from the dSph centre (analogous to rotation 
curves of spiral galax i es). Sp ecifically, we use the stellar kinematic 
data of I Walker et al.l ( |2009|) for the Carina, F ornax, Sculptor and 
Sextans dSphs, the data o f iMateo et al] ( l2008h for the Leo I dSph, 
and data from Mateo et al. (in preparation) for the Draco, Leo II 
and Ursa Minor dSphs. IWalker et al.l (I2009L W09 hereafter) have 
calculated velocity dispersion profiles from these same data un- 
der the assumption that l.o.s. velocity distributions are Gaussian. 



Here we re-calculate these profiles without adopting any particular 
form for the velocity distributions. Specifically, for a given dSph 
we divide the velocity sample into circular bins containing approx- 
imately equal numbers of member starsQ and within each bin we 
estimate the second velocity moment (squared velocity dispersion) 
as; 



N 



1 



{v)r 



(17) 



where A'^ is the number of member stars in the bin. We hold (V) 
fixed for all bins at the median velocity over the entire sample. 
For each bin we use a standard bootstrap re-sampling to estimate 
the associated error distribution for {V'^), which is approximately 
Gaussian. Fig. [9] displays the resulting velocity dispersion profiles, 
(V^)^^'^ (R), which are similar to previously published profiles. 

In order to relate these velocity dispersion profiles to dSph 
masses, we follow W09 in assuming that the data sample in each 
dSph a single, pressure-supported stellar population that is in dy- 
namical equilibrium and traces an underlying gravitational poten- 
tial dominated by dark matter. Implicit is the assumption that the 
orbital motions of stellar binary systems contribute negligibly to 
the measured velocity dispersions □ Furthermore, assuming spher- 
ical symmetry, the mass profile, A/(r), of the DM halo relates to 
(moments of) the stellar distribution function via the Jeans equa- 
tion: 



ld_ 

V dr 



(uv^) + 2 



/3(r)w2 _ GM{r) 



(18) 



where i'(r), Vr{r), and Pr = P{r) = 1 — Vg/Vr describe 
the 3-dimensional density, radial velocity dispersion, and orbital 
anisotropy, respectively, of the stellar component. Projecting along 
the l.o.s., the mass profile relates to observable profiles, the pro- 
jected stellar density I(R) and v elocity dispersion ap{R), accord- 
ing to l lBinnev & Tremainell2008l BT08 hereafter) 



(R) 



R^ 



i?2 



idr. 



(19) 



Notice that while we observe the projected velocity dispersion and 
stellar density profiles directly, the l.o.s. velocity dispersion pro- 
files provide no information about the anisotropy, /3(r). There- 
fore we require an assumption about /3(r); here we assume /3 = 
constant, allowing for nonzero anisotropy in the simplest way. 
For constant anisotropy , the Jeans equation has the solution (e.g., 
iMamon & LokaslbOOsI) : 



Gr 



^u{s)M{s)ds. 



(20) 



We shall adopt parametric models for I{R) and M (r) and then find 
values of the parameters of M (r) that, via Eqs. l |19t and l |20| l, best 
reproduce the observed velocity dispersion profiles. 



^ Kinematic samples are often contaminated by interlopers from the Milky 
Way foreground. F ollowing W09, we d iscard all stars for which the algo- 



rithm described by IWalker et al.l | |2009|) returns a membership probability 

le ss than 0.95. 

* lOlszewski et alj h996l) and iHargreaves et al.l l l 19961) conclude that this 
assumption is valid for the classical dSphs studied here, which have mea- 
sured velocity dispersions of ~ 10 km s^^. This conclusion does not 
necessarily apply to recently-discovered 'ultra-faint' Milky Way satel- 
lites, which have measured velocity dispersions as small as ~ 3 km s~^ 
iMcConnachie & C6tdl2010l) . 
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Figure 9. Velocity dispersion profile data for the 8 classical dSphs, obtained as described in the text (the impact of the binning choice is discussed in 
Appendix lHU . The solid lines correspond to the best-fit models for the inner slope when 7 is left free (dark), 7 is fixed to 1 (blue), and 7 is fixed to (red). 
Because of the large degeneracies among the halo parameters (see Section lsTl for a list), we do not list the corresponding best-fit parameters. The motivation 
for showing these profiles is to illustrate that our halo model is capable of describing the kinematic data, and that the inner profile is not constrained by the 
data. 



4.1.1 Stellar Density 



Stellar s urface densiti es of dSp hs are ty pically fit by IPlummej 
l l91lh. lKin3 h962l) and/ or ISersid mm . profiles (e.g., 



Ilrwin & Hatzidimitrioii|[l99^ . For simplicity, we adopt here the 



Plummer profile: 



[1 + R^/rl 



2 12 ' 

alfJ 



(21) 



which has just two free parameters: the total luminosity L and the 
projectetjl half-light radius rhaif- Given spherical symmetry, the 
Plummer profile implies a 3-dimensional stellar density (BT08) of: 



u{r) 



AI AR 



3L 



71- Jr di? /R2 



(22) 



Since we assume that DM dominates the gravitational potential at 
all radii (all measured dSphs have central mass-to-light ratios > 10, 
e.g., Mateo 19981), the value of L has no bearing on our analysis. 
We adopt values of rhaif (and associated errors) from Table 1 in the 
published errat um to W09; these data or i ginally come from the star 
count study of Ilrwin & Hatzidimitrioul ( Il995h . We have checked 



^ For consistency with I Walker et al ] 1 I2OO9I) we define rj^aif as the radius 
of the circle enclosing half of the dSph stellar light as seen in projection. 
Elsewhere this radius is commonly referred to as the 'effective radius'. 



that a steeper outer slope or a steeper inner slope for the light profile 
leaves unchanged the conclusions (see Appendix IH2l l. 



4.1.2 Dark matter halo 

For the DM halo we follow W09 in using a generalised Hemquist 
profile, as given by Eq. (|6). In terms of these parameters, i.e, the 
density ps at scale radius r^, plus the (outer, transition,inner) slopes 
{a, /3, 7), the mass profile is: 



M (r) = 47r 



s p{s)As = 

3-7 

'7 P-1 3 



3-7 



(23) 



- 7 + a 



where 2F1 (a, b; c; z) is Gauss' hypergeometric function. 

Eq. l|6j includes plausible halo shapes ranging from the 
constant-density 'cores' (7 = 0) that seem to describe rotation 
curves of spiral and low-surface-brightness galaxies (e.g., |de Blokl 
I2OIOI and references therein) to the centrally divergent 'cusps' 
(7 > 0) motivated by cosmological N-body simulations that model 
only the DM co mponent. For (q,/3,7) = (1,3, 1) E q. ([6j is just 
the cuspy NEW l lNavarro, Frenk & Whitell 199611 19971) profile. 
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4.2 Markov-Chain Monte Carlo Method 

For a given halo model we compare the projected (squared) velocity 
dispersion profile ap (R) (obtained from Eq. |19t to the empirical 
profile {V'^)(-R) (displayed in Fig.|9ll using the likelihood function 



c=n 



27rVar[(V2)(i?,)] 



: exp 



2 



Var[(l/2)(J?0] 



(24) 

where Var[{V'^){Ri)] is the variance associated with the empirical 
mean square velocity, as estimated from our bootstrap re-sampling. 

In order to explore the large parameter space efficiently, 
we employ Markov-chain Monte Carlo (MCMC) techniques. 
That is, we use the standard Met ropolis-Hastings algorithm 
etroDolis et al .ll953l : lHastingsll970h to generate posterior distri- 
butions according to the following prescription: 1) from the current 
location in parameter space, S„, draw a prospective new location, 
S', from a Gaussian probability density centred on Sn ; 2) evaluate 
the ratio of likelihoods at and 5"; and 3) if C{S')/C;{S„) > 
1, accept such that Sn+i = 5", else accept with probability 
C(S')/C(Sn) and S„+i = S„ with probability l-((S')/C{Sn). In 
order to account for the observati onal uncertainty associated wit h 
the half-light radius adopted from llrwin & Hatzidimitrio^ ( Il995l) . 
for each new point we scatter the adopted value of rhaif by a ran- 
dom deviate drawn from a Gaussian distribution with standard de- 
viation equal to the published error. This method effectively prop- 
agates the observational uncertainty associated with the half-light 
radius to the posterior distributions for our model parameters. 

Solutions of the Jeans equations are not guaranteed to corre- 
spond to physical models, as the associated p hase-space d i stribu- 
tion functions may not be everywhere positive. [An & Evansl ( l2006l) 
have derived a necessary relation between the asymptotic values of 
the logarithmic slope of the gravitational potential, the tracer den- 
sity distribution and the velocity anisotropy at small radii. Models 
which do not satisfy this relation will not give rise to physical dis- 
tribution functions. In terms of our parametrisation, this relation 
becomes 



7tr, 



2/?a 



(25) 



We therefore exclude from the Markov Chain those models which 
do not satisfy this condition. Because the Plummer profiles we use 
to describe dSph surface brightness profiles have 7traccr = 0, this 
restriction implies /3aniso J; 0. Given our assumption of constant 
velocity anisotropy, this disqualifies all radially anisotropic models. 
Relaxing this condition affects the results on the J-factors, but the 
difference is contained within their CLs (see Appendix IH2t . 

For this procedure we use the adaptive MCMC engine Cos- 
moMC jLewis & Bridld |2002|) . Although it was developed 
specifically for analysis of cosmic microwave background data, 
CosmoMC provides a generic sampler that continually updates the 
probability density according to the parameter covariances in order 
to optimise the acceptance rate. For each galaxy and parametrisa- 
tion we run four chains simultaneously, allowing each to proceed 
until the variances of parameter values across the four chains be- 
come less than 1% of the mean of the variances. Satisfying this 
convergence criterion typically requires ~ 10'' steps for our chains. 
We then estimate the posterior distribution in parameter space us- 
ing the last half of all accepted points (we discard the first half of 



available at |http://cosmologist.info/cosmomc| 



points, which we conservatively assume corresponds to the 'bum- 
in' period). 



5 DETECTABILITY OF MILKY WAY DSPHS 

This section provides our key results. For the benefit of readers who 
start reading here, we summarise our findings so far. 

In Section[2l we focused on generic (1, 3, 7) profiles, to show 
that, most of the time, the sub-structure contribution is negligible, 
and to check that the only relevant dSph halo parameters are the 
density nomalisation ps, the scale radius Va, and the inner slope 7 
(because JdSph oc rf^ x (Qintd)^^^"', see also AppendixlBt. 

In Section[3] we provided the sensitivity of present and future 
7-ray observatories, showing how it is degraded when considering 
'extended' sources (e.g. a flat profile for close dSphs), and an in- 
strument response that varies with energy. 

In Section|4l we presented our method to perform a Markov- 
Chain Monte Carlo analysis of the observed stellar kinematics in 
the 8 classical Milky Way dSphs under the assumptions of virial 
equilibrium, spherical symmetry, constant velocity anisotropy, and 
a Plummer light distribution. The analysis uses the observed ve- 
locity dispersion profiles of the dSphs to constrain their underlying 
dark matter halo potentials, parametrised using the five parameter 
models of Eq. (|6f. 



5.1 6-parameter MCMC analysis — varying 7 

Our kinematic models have six free parameters, for which we adopt 
unifoiTn priors over the following ranges: 

-logio[l-/3ani] : [-1,+!]; 
logio[ps/(Af0pc--^)] : [-10, +4]; 

Iogio[^=/Pc] : [0,4]; 

q: [0.5,3]; 
/3: [3,7]. 
7: [0,2] or [0,1]; 



The anisotropy parameter /3ani does not enter directly the pro- 
file/mass/J calculation, although it is of fundamental importance 
for the fit as it can correlate with the DM profile structure parame- 
ters (so with the mass and the J-factor). We have not checked ex- 
plicitly the details of these coixelations, but we have checked that 
restricting the range of possible /3ani does not significantly impact 
on the results for the J calculation. Hence, we do not discuss this 
parameter further below. 

5.1.1 Parameter correlations 

Fig. [To] shows the marginalised probability density functions 
(PDFs) of the profile parameters and the joint distributions of pairs 
of parameters. The features of these plots are driven by the fact 
that most of the stellar kinematic data lie at radii of up to few hun- 
dred parsecs (see Fig. [9). For instance, the outer slope j3 is not at 
all constrained (i.e. the fit is insensitive to the value of /?), because 
only tracers beyond a radius of r > 1 kpc are sensitive to this pa- 
rameter and these radii are sparsely sampled by the observations. 
The transition slope a and then the inner slope 7 are the two other 
least constrained parameters. In terms of best-fit models, as seen 
in Fig. [9] the match to kinematics data is equally good for varying 
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Figure 10. Joint distributions and marginalised PDFs of parameters entering the MCMC for the Draco dSph. The off-diagonal plots show joint distributions 
that highlight coiTelations between the parameters, while the on-diagonal plots are the marginalised PDFs of the parameters. This marginalisation includes the 
marginalisation over the velocity anisotropy parameter /3ani- (We do not plot a marginalised PDF or correlation for /3ani since it is a nuisance parameter for 
our analysis here.) 



7 (black) models and models in which we fix the value to 7 = 
(blue), or 7 = 1 (red). In the following, we will not discuss further 
the best-fit values. The more meaningful quantity, in the context of 
an MCMC analysis providing PDFs, is the median of the distribu- 
tion. 

Several groups have shown recently that in a Jeans anal- 
ysis, the observed flatness of dSph velocity dispersion profiles 
jWalker et al. leads to a constraint on A/(rhaif) — the mass 

enclosed within a sphere of radius rhaif — that is insensitive to as- 
sumptions ab out either anisotropy or the structura l parameters of 
the DM halo jWalker et al.ll2009l : IWolf et al.ll2010l) . Using for the 
appropriate radius the mass estimate Eq. (|9) and the above con- 
straint leads to a relation between the profile parameters 

log(ps) + 7log(rs) « constant. 

This relation explains the approximately linear correlations be- 
tween these parameters seen, for instance, in the bottom left panel 
of Fig.[Tl 



5.1.2 From p{r) to J(aint)-' uncertainty and impact o/7prior 

Fig. IFll shows the density profile for Draco as recovered by our 
MCMC analysis. It is noticeable that the confidence limits are nar- 
rower for radial scales of a few hundreds pc — this is a common 
feature of the density profile confidence limits for all the dSphs we 
have considered. As discussed above, this is partly due to the fact 
that these are the radii at which the majority of the kinematic data 
lie. The least constrained p{r) (less pronounced narrowing of the 
confidence limits) is that of Sextans, for which the range where 
useful data can be found is clearly the smallest compared to other 
dSphs (see Fig.|9). 

The variation of the constraints on p{r) as a function of ra- 
dius impacts directly on the behaviour of J. Complications arise 
because it is the profile squared that is now integrated along a l.o.s. 
(given the integration angle aint, see Eq.[5}- The median value and 
95% CL on J as a function of the integration angle Oint is plot- 
ted in Fig. [TT] (top), for two different priors on 7prior- The bot- 
tom panel gives the corresponding PDF for two integration angles. 
The prior has a strong impact on the result: the median (thick solid 
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Figure 11. Top: J(oint) for Draco as a function of the integration angle. 
Solid lines correspond to the median model, dotted lines to the 95% lower 
and upper CL. The to sets of curves correspond to two different 7pi ior f™ 
the MCMC analysis on the same data. Bottom: PDF of the J-factor for 
c'int = 0.01° (grey) and ajnt = 0.1° (black), when the range of the inner 
slope prior is [0 — 1] (solid lines) or [0 — 2] (dashed lines). 



curves and large symbols - top panel) is changed by ~ 50% for 
Oint <; 0.1°, but by a factor of ten for Uint ~ 0.01°. However, 
the most striking feature is the difference between the CLs: for the 
prior ^ 7prior ^ 2, the typical uncertainty is 3 to 4 orders of 
magnitude (red dotted curves), whereas it is only < than one order 
of magnitude for the prior < "fpr ior 1 (blue dotted curves) 
The bottom panel of Fig. [TT] shows that logjo J has a long and flat 
tail (associated to large 7 values). This tail is responsible for the 
large upper limit of the J-factor CLs for ^ 7prior 2. 

In Appendix |G2| a detailed analysis of the impact of these 
two priors is carried on artificial data (for which the true profile 
is known). We find that the prior ^ 7prior ^ 2 satisfactorily 
reconstructs p(r) and J(aint), i.e. the MCMC CLs bracket the true 
value. This is also the case when using the prior ^5 7prior ^ 1- 
However, two important points are noteworthy: 



Note that this behaviour is grossly representative of all dSphs, although 
the integration angle for which the uncertainty is the smallest, and the am- 
plitude of this uncertainty depend, respectively, on the dSph distance (see 
Section |23] for the generic dependence), and on the range/precision of the 
kinematic data (see above). 



• this prior obviously performs better for ^ 7truc ^ 1 profiles 
where it gives much tighter constraints on J; 

• for cuspier profiles (e.g., 7truc = 1-5), this prior succeeds 
slightly less (than the prior ^ 7prior ^ 2) in reconstructing 
p(r), but it does surprisingly better on J in terms of providing a 
value closer to the true one (see details and explanations in Ap- 
pendix [G2](- 

DM simulations and observations do not favour 7 > 1, al- 
though steeper profiles can still fit the kinematic data in a Jeans 
analysis. Indeed, the Aquarius simulations indicate values of 7 
slightly smaller than 1, and although some recent simulations 
jlshivama et al.|[201ol) have argued for cuspy profiles, this happens 
for micro-haloes only. Given that the J-factor for the cuspier pro- 
files are only marginally more (or even less) reliable when using the 
prior ^ 7prior ^ 2, we restrict ourselves to the ^ 7prior ^ 1 
prior below. 

Note that other sources of bias exist. First, the reconstruc- 
tion of p{r) or J(aint) is affected by the choice of binning used 
in the estimation of the empirical velocity dispersion profiles. Ap- 
pendix |HT] shows that we obtain slightly different results when we 
apply our method to empirical velocity dispersion profiles calcu- 
lated from the same raw kinematic data, but using different num- 
bers of bins. We find that the effects of binning add an extra factor 
of a few uncertainty on J for the least well measured (in terms of 
radial coverage) dSphs, for which more measurements are desir- 
able. (On the other hand, Fornax and Sculptor are found to provide 
robust results against different binnings.) Second, we note that the 
analysis presented here uses a fixed profile for the light distribution 
which, when combined with our assumption of constant velocity 
anisotropy, restricts the possible halo profiles we can recover. Our 
constraints on p{r) a nd J(aint) are there fore sensitive to these as- 
sumptions (see, e.g., IStrigarietal.ll2010l for an example of fitting 
the dSph kinematic data with cusped profiles when the light pro- 
file is also allowed to be cusped), although this does not change 
our conclusions (see Appendix IH2I where different light profiles 
are used). This situation is set to change over the coming years as 
new distribution function-based models will permit constraints to 
be placed on the slope of the DM density profiles (Wilkinson et al. 
201 1, in prep.). 



5.7..? Best constraints on J: median value and CLs 

As validated by the simulated data, we are now able to provide ro- 
bust (although possibly not the best achievable with current data) 
and model-independent constraints on J(aint) for the 8 classi- 
cal dSphs. The results are summarised in Table |2] in terms of the 
median, and 68% and 95% CLs. The J-factor is calculated for 
Q!int = 0.01° (an angle slightly better than what can be achieved 
with FCA), Qint = 0.1° (typical of the angular resolution of exist- 
ing GeV an d TeV 7-ray instru ments), and for Uc — 2rhaif/rf (as 
proposed in I Walker et al.ll2oTll) . We do not report the values of ps 
and Vs as these vary across a large range — and therefore do not give 
additional useful information — nor the value of 7 as it is forced in 
the range ^ 7prior ^ 1 to give the least biased J value. 

There is no simple way to provide unambiguously the best 
target, as their relative merit depends non trivially on their dis- 
ta nce, their ma s s and the integration angle selected. As proposed 
in lWalker et al. ll), since the most robust constraint on J is ob- 
tained for Oint = Qc, having different integration angles for each 
dSph can be a good starting point to establish a relative ranking. 
The situation is complicated further for background-limited instru- 
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Table 2. Positions of the classical dSphs iMated Il998t) sorted according to their distance: longitude, latitude, distan ce, 2rv,3if (taken from 
llrwin & Hatzidimitriou|[l995l) , the galactic angle away from centre 4> = cos~^ [cos(long.) cos(lat.)], and ~ 2rYi^\f/d (see lWalker et aljboi ll) . The 
remaining columns are the median value with 68% (95%) CLs for A/300 ™d logjg [ J(aint)] from the six-parameter MCMC analysis (0 ^ 7prior ^5 !)■ For 
conversion factors to units used in other studies, please refer to numbers given in Appendix IaI 
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* Note that the values for logj^Q [J(cic)] differ from those quoted in lWalker et al.l ( l201 lb as the MCMC analysis is slightly different here. 



merits such as CTA, as some lo ss of sensitivity can occur (see, e.g. 
Figure 4 of IWalker et 3111201 l|). This is discussed, taking into ac- 
count the full detail of the instruments, in Section 15.31 However, 
in this respect, the best target for future instrument may eventually 
become Leo II, which despite a quite large uncertainty outshines 
all other dSphs at aint = 0.01° (see also Fig.ll2t. We note how- 
ever that it is the dsph with the smallest amount of kinematic data 
at present (so it has the most uncertain J-factor). 



5.1.4 dSphs in the diffuse galactic DM signal: contrast 

The uncertainties in J are illustrated from a different viewpoint in 
FigdU It shows, in addition to the mean, 68% and 98% CLs on 
the Js, the latitudinal dependence of the Galactic DM background 
(smooth and galactic clump contribution) for the same integration 
angle[3 For a typical present-day instrument resolution (integra- 
tion angle aint ~ 0.1°), we recover the standard result that the 
Galactic Centre outshines all dSphs. 

The three panels illustrate the loss of contrast (signal from the 
dSph w.rt. to the diffuse Galactic DIM signal) as the integration 
angle is increased. This is understood as follows: the integrand ap- 
pearing in Eqs. iC4) and iC5\ is mostly insensitive to the l.o.s. di- 
rection a few tens of degree away from the Galactic centre, so that 
Eq. l lC6t holds, giving an af^t dependence. 

For detectability (see also Sec[3j, the naive approach of max- 
imising the integration angle (to maximise JdSph) must be weighed 
against the fact that an increased integration angle means more 



The smooth profile is taken to be an Einasto profile, the clump dis- 
tribution is a core one, where as their inner profile are Einasto with con- 
centration and parameters a /a lBuUock et alj j200lh Normalising the mass 
distribution to have 100 clumps more massive than lO*Af0, and taking 
dP/dM QC M~'^-^ leads to a DM fraction into clumps o f ~ 10% for 
clumps distributed in the range 10"'^ - 10^°Mq (see, e.g.. lLavalle eTal] 
I2OO8I and references therein). The local DM distribution is fixed to the fidu- 
cial value pQ = 0.3 GeV cm^'^. The exact configuration is unimportant 
here as this plot is mostly used for illustration purpose. 



astrophysical 7-ray and cosmic-ray background. For large inte- 
gration angles, dSphs also have poor contrast against the diffuse 
Galactic DM annihilation signal, indicating that the Galactic halo 
is a better target for any search on angular scales >I (see e.g. 



my 

lAbramowski et al.l201Ibl for such a search with H.E.S.S.) 



5.1.5 Comparison to other works 

Comparison between different works can be difficult as every au- 
thor uses different definition, notations and units for the astrophys- 
ical factor To ease the comparison, we provide in Appendix lAlcon- 
version factors between standard units (we also point out issues to 
be aware of when performing such comparisons). 

Below is a comparison to just a few of the works published 
on the subject, and only for the objects that these studies and the 
present one have in common: 

• The lEvans et alj l l2004h values of J/ASl for Draco (with 
An — 10~^ i.e. Clint = 0.1°) for all the profiles they explored 
(cored, 7 — 0.5, 7 = 1, 7 = 1.5) are larger (after correction by 
AQ, given their definition of the astrophysical factor) than our 95% 
CL upper limit for this object shown in Table |2] The difference is 
probabl y related to our data set which is about twice as large as that 
used bvlsvans et al lll20()4l) . 

• IStrigari et al. "(12007 ) provide directly the 7-ray flux (i.e. in- 
cluding the particle physics term), so that we can only compare our 
respective rankings. These agree in general but for Sculptor we find 
a lar ger flux than Draco, conversely to these authors. 

• IPieri e t al. (2009) focused on Sextans, Carina, Draco and Ursa 
Minor. They found the latter to have the largest J ("I>cosmo in their 
notation) of these 4 objects, followed by Draco, Carina and Sex- 
tans. But for the last two, this ranking is similar to ours. However, 
while their values of J fall within our 68% (UMi, Sextans) or 95% 
(Car ina) CL, the i r value for Draco is above our 95% CL upper limit. 

• lEssig et all ( I2OO9I) also performed a statistical study on Draco 
and Ursa Minor, to determine their profiles from kinematic data 
and to derive the confidence levels on the J-factor. Given that their 
integration is performed on a slightly larger opening angle (0.14°), 
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Figure 12. Galactic contributions to J for the smooth (blue-dashed line), 
mean clump (red-dotted line) and sum (black-solid line) vs the angle from 
the Galactic centre. The symbols show J for the dSphs, assuming a prior of 
^ 7prior ^ 1 on the central DIVI slope. The central point con'esponds to 
the median values, the solid bars to the 68% CLs, and the dotted bars to the 
95% CLs. The integration angle is, from top to bottom 0.01°, 0.1°, and 1°. 
The Galactic contributions Jsm and ( Jsubcl) scale as a^-jj, but JdSph does 
not, changing the contrast of the dSphs w.rt. to the DM Galactic background 
(see text for details). 



compared to the median and confiden ce levels we de rived in Ta- 
ble [2] for aint = 0.1°. The values of lKuhlenl ilOld) (multiplied 
by 4n to match our definition of J) generally fall inside our 68% 
CL intervals, but for Leo II his value is just within our 95% CI 
while Draco and Carina c annot be accommodated at all. For these 
two objects, the values o f iKuhlerj J201ol) are much larger than the 
ones we find, and this is unlikely to be explained by the 0. 05° dif- 
ference in integration angles. A simple explanation is that iKuhleij 
C2010i') does not use stellar kinematic data directly in his analysis, 
but stacks suitable Via Lactea halos (M300 ~ IO^Mq and appro- 
priate distances) and uses those averages to estimate J. Focusing 
on the ranking (without worrying about contrast to the background 
and the other instrumental constraints), both we and lKuhlenlfcoid) 
agree that among the classical dSph UMi is a most promising tar- 
get. However, wh ile we find Scu lptor and Draco to be the next most 
favorable targets . iKuhlenI JlO 1 Ot) names Draco and Carina from his 
'simulation-based' approach. 

For completness, we also compare our median values with the 
J values used by different experimental groups: 

• The IVIAGIC collaboration published point source lim- 
its for Draco d Albert et a l. 2008) adopting the scheme of 



ISanchez-Conde et al. TiootI) of a power law density profile, with 
an exponential cut-off. They examine two scenarios, a cored and a 
cusped model, but find no discernable difference when calculating 
J for integration angles < 0.4°, i.e. larger than the IVIAGIC PSF. 
The value of J they calculate for Draco is a higher than ours (after 
appropriate scaling of the integration region and unit conversion) 
by about a factor of 2. 

• The VERIT AS collaboration als o published limits on Draco 
and Ursa IVIinor jAcciari et al.ll2010|). They assume a NFW pro- 
file, take the density pr ofiles from Strigari et al.l (l2007b and follow 
ISergstrom et al.l ( Il998h for the calcula tion of J. Whilst the range of 
density values in 
the values used in 



Strigari et al 



Acciari et al 



( 20071) have a physical motivation 
(20103 are rather arbitrarily chosen 
to be the midpoint of that range, which leads to consistently higher 
J- values than ours (by a factor of 3 for Draco and a factor of 1.2 
for UMi) . 

• The H.E.S.S. collaboration jAbramowski et aLll201 lal) pub- 
lished limits on the southern sources Sculptor and Carina using 
NFW and isothermal profiles with a number of varying assump- 
tions. This leads to a range of calculated J-values (rather than a 
single solution) that are consistent with our median value and esti- 
mated uncertainties. 

• The Fermi collaboration jAbdo et al .I2OIOI) has published lim- 
its for a number of the sources studied here. They adopted a NFW 
profile within the tidal radius and following llVIartinez et al.l l [2009l) 
they calculated the J- value (using an IVICMC approach on the ob- 
served stellar velocities) for a 1° integration angle which is compat- 
ible with their high energy PSF. From this they find Draco to have a 
larger J compared to the other dwarfs (a factor of ~ 2 higher than 
the next dwarf which is Ursa IVIinor), contrary to what we find in 
this study. 



our results appear to be in agreement. Their 90% CL limits are 2-3 
times larger than the 95% CL limits given in Table |2]but this may 
be due to the larger range they adopt for the prior on the inner slope 
(see App.|Gl| . 

• IKuhlerj (201(3) gives the astrophysical factors of all the dSphs 
using a point-like approximation and a NFW DIVI profile, and in- 
tegrated with a Qint = 0.15° angular resolution. These can be 



5.2 5-parameter MCMC analysis: 7prior fixed 

Higher resolution numerical simulations following both DM 
and gas, additional kinematic data and new modelling tech- 
niques may help constraining the value of 7 in the near fu- 
ture. With the knowledge of 7, we should better constrain the 
radial-dependence of J, which is crucia l to disentangle, e.g . 
dark matter annihilation from DM decay dBovarskv et al.ll2006l : 
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omares-Ruiz & Siegal- GaskinsI lloToh . The topic of decaying 
dark matter goes beyond the scope of this paper, and it will be dis- 
cussed elsewhere. Below, we merely inspect the gain obtained on 
the J prediction when having a strong prior on 7, and briefly com- 
ment on the possibility to disentangle 7 = profiles from 7 = 1.0 
profile in the case of annihilation (if this cannot be achieved, hopes 
for disentangling decay from annihilation would be quite low on a 
single object). 

5.2.1 Parameter correlations 

We repeat the MCMC analysis for fixed value of the inner slope 
7prior = 0., 0.5, 1., and 1.5. The priors for the five other parameters 
are as given in Section lSTI 

Usin g Eg. lO for the mass having a robust estimate of 
Mfr-haif ) Iwalker et al.l2009l : IWolf et al.l201ol : lAmorisco & EvansI 
l201lh gives log(pB) + 7log(ra) « constant which reduces to 
log(ps) ~ constant for 7 = 0. As a result, we expect a strong 
correlation between ps and Tb when 7prior = 1 and none when 
7prior = 0. This Is Confirmed by the result of our MCMC analysis 
shown in Fig. [T3] (here, for the Draco case). The half-light radius 
Thaif for Draco is ~ 200 pc, but we choose to show the PDF for 
M300 in the bottom panel of Fig. [13] as we wish to compare the 
mass of the dSphs among themselves (see Table|2j. It confirms that 
the mass within an appropriate radius can be reliably constrained 
by the data regardless of the value of 7. 

5.2.2 Uncertainties on the profile and on J 

For any given 7, the uncertainty on p(r) at small radii is related to 
the range of Tb values at which the asymptotic slope is reached (for 
each profile accepted by the MCMC analysis). For 7prior = 0, the 
maximum uncertainty on p{r) is directly related to the maximum 
uncertainty on ps (since for r <Si rs, p(r) is constant) which can 
be read off the PDF (top-left panel of Fig. [T3l l. This leads to an 
order of magnitude uncertainty on p{r) for small r, which is con- 
sistent with the 95% CL shown in top panel of Fig. [14] For 7 > 0, 
the uncertainty has to be read from the dispersion in the values of 
PsrJ, or equivalently, the mass Mzw- The bottom panel of Fig. 1131 
shows that this mass is well-constrained, independently of 7 for the 
case of Draco (see however in Table|2]for a larger spread for some 
dSphs), resulting in a smaller uncertainty for 7prior = 1-5 than for 
7prior = (top panel of Fig.[T3ll. We checked that the CLs obtained 
in Fig.[T4](in Appendix|G2) for the artificial data enclose correctly 
the range of reconstructed values: they are consistent with a larger 
reconstruction bias for 7prior = than for 7prior = 1.5 at small 
radii. 

For the uncertainty on J, we can obtain a crude estimate by 
relying on the approximate formulae given in Appendix |B] For 
7 > 0, J oc plrl, and substituting the constant A/goo relation- 
ship leads to J oc rg"^"'. The value of Ts, as seen in its PDF in 
the top and middle panels of Fig.[T3] varies by roughly a factor of 
10. Because of the weighting power 3 — 27, the uncertainty on J is 
expected to be the smallest for 7 = 1.5, which is in agreement with 
the curves in Fig. [14] (bottom panel). However, the analysis of the 
artificial data in Appendix IG2 1 shows that the typical CL on J ob- 
tained in the bottom panel of Fig.[T4]is likely to be underestimated 
for 7prior = 1-5 (up to factor 0(2), see Fig. [G2J-E3 This happens 
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Figure 13. Top: correlation and PDF of the profile parameters ps and rs 
from the 5-pai'ameter MCMC analysis 7pi-ior = 0.0. Middle: same, but 
for 7prior = l-O- Bottom: PDF of M300, the mass at 300 pc. 



for any integration angle. For this reason, we cannot rely of the J 
value for 7prior = 1.5 and focus only on the three cases 7prior = 0, 
7prior = 0.5, and 7prior = 1.0 below. 



" This is understood as for the latter, the inner region (r <C ) contribute 
the most to J, and even small differences for p(r ~ rs) are bound to trans- 



late in sizeable differences for p(r — > 0). Conversely, similar differences on 
p for shallower profiles is not an issue as their inner parts do not contribute 
to J. 
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Figure 15. Median J-factor values (symbols) and 68%/95% CLs (solid bars; dashed bars) for the fixed 7prior analysis (the result for 7pi.ior = 1.5 is not 
shown because it is not reliable, see Sect. lG2t . The blue dashed fine shows the expected scaling with distance for point sources: 3.1 • 10^^ ci~'^ I^q kpc~^]. 
The panels show, from top to bottom, three integration angles aint = 0.01°, 0.1°, and Oc ~ Sr^/d (an angle very similar to the angle enclosing 80% of 
the flux, see Fig. ll6t that optimises the detennination of the J-factor for a given dSph (hence the error bars are smaller in this plot than in the other two). The 
yellow solid fines (and broken lines in tfie bottom panel) correspond to tfie Galactic DM background including botfi tfie smootfi and clumpy distiibutions. For 
the bottom panel, this is not a smooth curve since it depends on the integration angle ajnt that varies from dSph to dSph in this figure. Note that the choice of 
using the critical angle ai^t = etc is optimal in the sense that it gives the most constrained value for J. But where the Galactic background annihilation signal 
approaches that of the dSphs (see for example. Sextans and Fornax), the motivation for staring at the dSphs rather than simply looking at the Galactic halo is 
gone. 



5.2.3 J{d) and departure from the 1/d? scaling 

Fig-EDshows the J median values, 65% and 95% CIs as symbols, 
daslied and solid error bars respectively, for an inte gration angle 
of 0.0 1° (top), 0.1° (middle), and Qc « 2rhaif/d ( I Walker et all 
1201 ih The X-axis is the distance to the dSph (in kpc). For point-like 
sources, the J-factor of a single dSph scales as as illustrated 
by the blue-dashed line. Departure from this scaling is interpreted 
as a combination of a mass effect and/or a profile effect. For in- 
stance. Sextans and Carina are dSphs with smaller M300 with re- 
spect to the other ones (see Tab. |2j; consequently they are located 
below the dashed blue line in the top panel of Fig. [15] The excep- 
tion is Leo II, which has a 'small' mass but is nevertheless above 



the dashed line. Although this analysis cannot constrain 7, we are 
tempted to interpret this oddity in terms of a 'cuspier' profile (w.r.t. 
those for other dSphs), which would be consistent with the fact that 
its J remains similar in moving from aint = 0.1° (middle panel) 
to 0.01° (top panel). However, an alternative explanation (which 
would be more consistent with the results obtained in this paper) 
could be the fact that Leo II has the smallest amount of kinematic 
data at present, and that its J is overestimated (see Appendix IHll to 
support this line of argument). We repeat that the relative brightness 
of the dSphs is further affected for background-dominated instru- 
ments (as described in Sec. |3}, so that the ranldng has to be based 
on Fig[T6]discussed in the next section. 
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Figure 14. Median values (solid lines, filled symbols) and 95% CLs (dashed 
lines, empty symbols) from the fixed 7prior MCMC analysis on Draco. 
Top: density profiles (the gray arrow indicates the value of r^aif ). Bottom: 
J-factor (the gray arrow indicates etc ~ 2rhaif/rf). 



The bottom panel of Fig. [15] shows the J value for an 'opti- 
mal' integration angle Oc that is twice the half-light radius divided 
by the dSph distanc43 (this corr esponds to the inte gration angle 
that minimises the CLs on J; see lWalker et al.ll201 ll) . The yellow 
broken solid lines show the expected signal from the diffuse Galac- 
tic DM annihilation background, including a contribution from 
clumpy sub- structures (the extragalactic background, which also 
scales as apnt^ has not been included). The total background may 
be uncertain by a factor of a few (depending on the exact Galactic 
(smooth) profile and local DM density). Its exact level — which de- 
pends on the chosen integration angle — determines the condition 
for the loss of contrast of the dSph signal, i.e. the condition for 
which looking at the DM halo (rather than at dSphs) becomes a 
better strategy. 



5.2.4 Conclusion for the fixed 7prior analysis 

The analysis of simulated data shows that the analysis for 7piior ~ 
1.5 is biased by a factor of O{10) and that the CLs obtained on the 
real data are likely to be severely under-estimated in that case. But 
such steeply cusped profiles are neither supported by observations 
nor motivated by current cosmological simulations. For values of 



CLs for J(aint ) are provided along with the paper for readers interested 
in applying our analysis to existing and future observatories. 



7prior ^ 1, this bias is a factor of a few only, so that it shows that 
the results from a fixed 7prior analysis of the 8 classical dSphs are 
robusts. However, this analysis shows that unless very small inte- 
gration angles aint < 0.01° are chosen (or if 7truc 1), knowing 
the exact value of 7 does not help in improving the determination 
of J. Indeed, even using Draco, the stellar population of which is 
one of the most studied, the CLs of the three reconstructed fluxes 
(7piior = in black full circles, 7prior ~ 0.5 in red triangles, and 
7prior = 1-0 in blue stars) in Fig.[T4](bottom), overlap. Reversing 
the argument, if we do not know the inner slope, and if a 7-ray sig- 
nal is detected from just one dSph in future, there will be little hope 
of recovering the slope of the DM halo from that measurement only. 

This means that the best way to improve the prediction of the 
J-factor in the future relies on obtaining more data and a more 
refined MCMC analysis; an improved prior on the DM distribution 
makes little difference. 



5.3 Sensitivity of 7-ray observatories to DM annihilation in 
the dSphs 

The potential for using the classical dSph to place constraints on 
the DM annihilation cross-section, given the uncertainties in the 
astrophysical J-factor, can be seen in Fig. [T6] Previous analyses 
have adopted the solid angle for calculation of the J-factor to be 
the angular resolution of the telesc ope for a point-like source, typi- 
cally assuming a NFW-like profi le jAcciari et al.ll2010l : lAbdo et al.l 
I2OIOI : lAbramowski et al]|201 lal) . By contrast our sensitivity plots 
take into account finite size effects: i) the J values are based on the 
MCMC analysis with the prior ^ 7prior ^ 1, where the corre- 
sponding J are shown in Fig. [12] ii) the energy dependent angular 
resolution has also been taken into account assuming a standard 7- 
ray annihilation spectrum (see Section |2.1.U . Moreover for Fermi- 
LAT the background level assumed has been increased (resulting 
in a 25% worsening of the sensitivity above 100 MeV) to reflect 
the average situation in the directions of the classical dSph (the 
variation between the individual dSph is only 7% rms). A likeli- 
hood based analysis is used for both FCA and Fermi and a nominal 
observation zenith ang le of 20° assumed for the FC^EI 

(see Sec- 

tionim. 

The panels from top to bottom coiTespond to increasing DM 
(neutralino) masses. At low values, Fermi has a better sensitivity 
than FCA; at a mass of about 1 TeV the two are comparable, and 
for higher masses the FCA becomes the more sensitive instrument 
due to the vastly greater effective area at the photon energies at 
which the annihilation spectrum is expected to peak. Note that the 
precise value of (av) where the relative sensitivities of the two in- 
struments cross depends on the form of the DM annihilation spec- 
trum. Since we are examining the uncertainties in the astrophysical 
J-factor to the detectability of dSphs, we have used a conservative 
spectrum averaged over a number of possible annihilation channels 
(see Fig.[T]( which results in the majority of produced 7-ray pho- 
tons having energies ~10% of the DM particle mass. If we were 
to move from a relatively soft spectrum, such as bb to a harder 
one, such as t^t^ , this would benefit both instruments in differ- 
ent ways. For Fermi-LAT a harder spectrum makes the signal eas- 
ier to distinguish above the diffuse 7-ray background; indeed the 



The energy threshold for a ground based instrament is dependent on the 
zenith angle of observation. This means that the actual energy threshold for 
a given object will depend on the object's declination and the latitude of the, 
yet to be determined FCA site. 
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Figure 16. Minimum detectable < av > for known dSphs shown as a function of their distance, for different assumed DM masses (separate panels). 100 hour 
observations with an FCA (red circles) are compared to 5 years of Fermi observations (blue squares). Error bars indicate 65% (solid lines) and 95% (dotted 
lines) confidence limits. The integration angle is adapted to the Qc of each dSph and the energy-dependent PSF of the two instruments. The strategy of using 



ago, rather than Oc, is indicated with hollow symbols for the FCA case. The line for < av >= 
the panels. 



10- 



' cm^ s ^ is drawn for compaiison purpose between 



lAbdo et all JlOlOh found that the detectable flux limit from a po- 
tential source could vary by a factor of 2-20 (with lower particle 
masses benefiting the most) between these different annihilation 
spectra. For the FCA, which has a very large effective area to pho- 
tons ^ 100 GeV, the benefits of having more high energy photons 
is very apparent when it comes to flux sensitivity. For both obser- 
vatories, an increased number of high energy photons needs to be 
balanced with the correspondingly better angular resolution, partic- 
ularly if (e.g. for Fermi-LAT) a point-like source becomes spatially 
resolved. 

Our analysis places Ursa Minor as the best candidate for the 
northern sky (marginally better than Draco, which has long been a 
favourite target of northern hemisphere observatories) and Sculp- 
tor for the southern sky, when it comes to a favourable median and 
low uncertainty in the J-factor. It should be noted, however, that 
although the closest objects seem to be favoured, Leo II has the po- 



tential to yield a stronger signal, however more kinematic data are 
needed in order to constrain better its I-factor. In addition, it should 
be noted that the uneven sensitivity of the Fermi-LAT across the 
sky, caused in particular by the proximity of bright source^3 as 
well as the galactic diffuse background can change what is consid- 
ered the favorite candidate. 

We emphasise that in our analysis the inner slope 7 has not 
been constrained, but that a better independent determination of 
7 in future will not help providing a better determination of J (see 
Fig. list; this is discussed further in the Appendices. Carina, Fornax 
and Leo I are the targets least favoured. W hen compared to existing 
limits from Fermi-LAT jAbdo et^ .I2OIOI) or the current generation 

1^ In particular there is a bright GeV emitter IFGL J0058.4-3235 only 
~ 1.1° away from Sculpt or which significan tly worsens the upper limit on 
that object as discussed bv lAbdo et aljj201ol) . 
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Figure 17. Sensitivity reach in tlie — (crv) plan for FCA (100 hr) and Fermi (5 yr), for our best candidates UMi (median value and 95% CIs) and Leo II 
(median only). Black asterisks repr esen t points from MSSM models that fall within 3 standard deviations of the relic density measured in the 3 year WMAP 
data set (taken from lAcciari et al..2010l) . 



of ACTs jAcciari et allboiol : lAbramowski et al.llZOllal) it can be 
seen that our limits are not dissimilar from those that have already 
been published. For Fermi this is not surprising, since the source 
is unresolved and any difference should relate only to the assumed 
increase in exposure from 1 to 5 years, resulting in a factor of a 
few at best. The similarity in sensitivity between current and future 
ACTs is perhaps more surprising, but this as stated earlier relates 
to the naive assumptions made on the form for the J-factor and the 
solid angle integrated over; in order to reach the currently claimed 
limits requires a deep exposure with an instrument as sensitive as 
CTA. 

One last thing to note is that a common way to synthesise a 
deeper exposure is to stack observations of different sources to- 
gether to provide an effective long exposure of a generic source. 
For a common universal halo profile this may be fine, however any 
analysis will have to take into account the different integration an- 
gles for each individual source correctly. If all dSphs do not share 
a common halo profile and hence have different 7 values, we have 
to rely on the varying-7 analysis presented in the previous section 
and the relative ranking of potential targets would then be different. 



6 DISCUSSION AND CONCLUSIONS 

We have revisited the expected DM annihilation signal from dSph 
galaxies for current (Fermi-LAT) and future (e.g. CTA) 7-ray ob- 
servatories. The main innovative features of our analysis are that: 
(i) We have considered the effect of the angular size of the dSphs 
for the first time. This is important since, while nearby dSphs have 
higher 7 ray flux, their larger angular extent can make them sub- 
prime targets if the sensitivity is limited by cosmic ray and 7-ray 
backgrounds, (ii) We determined the astrophysical J-factor for the 
classical dSphs directly from photometric and kinematic data. We 
assumed very little about their underlying DM distribution, mod- 
elling the dSph DM profile as a smooth split-power law, both with 



and without DM sub-clumps, (iii) We used a MCMC technique to 
marginalise over unknown parameters and determine the sensitiv- 
ity of our derived J-factors to both model and measurement un- 
certainties, (iv) We used simulated DM profiles to demonstrate that 
our J-factor determinations recover the correct solution within our 
quoted uncertainties. 

Our key findings are as follows: 

(i) Sub-clumps in the dSphs do not usefully boost the signal. 
For all configurations where the sub-clump distribution follows the 
underlying smooth DM halo, the boost factor is at most ~ 2 — 3. 
Moreover, to obtain even this mild boost, one has to integrate the 
signal over the whole angular extent of the dSph. This is unlikely 
to be an effective strategy as the diffuse Galactic DM signal will 
dominate for integration angles aint 1°- 

(ii) Point-like emission from a dSph is a very poor approxi- 
mation for high angular resolution instruments, such as the next- 
generation CTA. For a nearby dSph, using the point-like approxi- 
mation can lead to an order of magnitude overestimate of the de- 
tection sensitivity. In the case of a nearby cored profile consisting 
of very high mass DM particles, a point source approximation can 
be unsatisfactory even for the modest angular resolution of Fermi- 
LAT. 

(iii) With the Jeans' analysis, no DM profile can be ruled out by 
current data. The use of the MCMC technique on artificial data also 
shows that such an analysis is unable to provide reliable values for 
J if the profiles are cuspy (7 = 1.5). However, using a prior on the 
inner DM cusp slope ^ 7prior ^5 1 provides J-factor estimates 
accurate to a factor of a few. 

(iv) The best dSph targets are not simply those closest to us, as 
might naively be expected. A good candidate has to combine high 
mass, close proximity, small angular size (< 1°; i.e. not too close); 
and a well-constrained DM profile. With these criteria in mind, we 
find three categories: well-constrained and promising (Ursa Minor, 
Sculptor and Draco), well-constrained but less promising (Carina, 
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Fornax and Leo I), and poorly constrained (Sextans and Leo II). 
Leo II may yet prove to be a viable target as it has a larger median 
J-factor than UMi, however more data are required to confirm its 
status. 

(v) A search based on a known DM candidate (from, e.g., forth- 
coming discoveries at the LHC) will do much to optimise the search 
strategy and, ultimately, the detection sensitivity for all 7-ray obser- 
vatories. This is because the shape of the annihilation spectrum is 
a strong driver of the photon energy range that can provide the best 
information on the candidate DM particle mass. Fermi-LAT has 
great potential to probe down to the expected annihilation cross- 
section for particles of mass 700 GeV, whereas a ground based 
instrument is more suited for probing particle masses above a few 
hundred GeV with a sufficiently deep exposure. However, even for 
5 yr of observation with Fermi-LAT or 100 hrs with FCA, the sen- 
sitivity reach (Fig. I17t remains anywhere between 4 to 10 orders 
of magnitude above the expected annihilation cross-section for a 
cosmological relic (depending on the mass of the DM particle can- 
didate). Improving these limits will require a harder annihilation 
spectrum than the conservative average we have adopted in this 
study, or a significant boost (e.g. from the Sommerfeld enhance- 
ment) to the 7-ray production. 

Finally, the ultra-faint dSphs have received a lot of interest in 
the community lately, as they could be the most-DM dominated 
systems in the Galaxy. We emphasise that the MCMC analysis we 
have performed for the classical dSphs cannot be applied 'as is' for 
these objects. First, the sample of stars observed is smaller. Sec- 
ond, the velocity dispersion is smaller and suffers from larger un- 
certainties than those for the classical dSphs. The robustness and 
systematic biases of the MCMC analysis will be discussed else- 
where (Walker et al., 2011, in preparation). Results concerning J 
for the ultra-faint dSphs will be presented in a companion paper. 
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APPENDIX A: DEFINITIONS, NOTATION, CONVERSION 
FACTORS 

Studies of DM annihilations in the context of dSphs involves both 
particle physics and astrophysics. The obvious difference of scales 
between the two fields and habits among the two communities have 
given rise to a plethora of notations and unit choices throughout the 



literature. In this Appendix, we provide some explanatory elements 
and conversion factors to ease comparison between the different 
works published on the subject. 

As mentioned in ij2l we define the differential 7— ray flux as 
integrated over the solid angle Afi as 

^(S^,An)=<I>PP(iJ^)x J(Af7), 

where 

$pp(£; ) = 5!£i = _Li£±fifi!i £^ 

and 

J(An) = / [ plyi{l,^)dldQ. 
J An J 

The solid angle is simply related to the integration angle Qint by 

AQ = 27r • (1 - cos(Qint)) • 
In our work, the units of these quantities are as follows: 

• [d$^/d£;^] = s"^ GeV"^ 

• [c&PPCE'^)] = cm^ GeV-^( sr"^); 

• [J] =Mgkpc-5(sr). 

First of all, note that the loca tion of the l/47r f actor appearing in 
t&PP is arbitrary. We followed [Pieri et al.l (|2009|) and included it in 
the particle physics fac tor. In other works, it ca n appear in the astro- 
physical factor J (e.g.. lBringmann et al.ll2009l) . Therefore, to com- 
pare the astrophysical factors between several studies, one must 
first ensure to correct the value of J by 47r if needed. In the text, 
we did not explicitly stated the solid angle dependence in the units 
of J as it is dimensionless quantity. The conversion factor (once 
the 47r issue is resolved) from our J units to that traditionally found 
in the literature are: 

• 1 Mq kpc-^ = 10"" A4 

• 1 m| kpc"-"^ = 4.45 X 10'' GeV^ cm^^ 

• 1 kpc"^ (sr) = 1.44 x lO"^"^ GeV^ cm"'* kpc (sr) 

Before comparing any number, one must also ensure that the solid 
angle ASl over which the integration is performed is the same. In 
most works, a Qint ~ 0.1° angular resolution is chosen, corre- 
sponding to AO, — 10^^ sr. However this is not always the case, as 
in the present study where we explore several angular resolutions. 
Note that the quantity J = J/ASl (in GeV^ cm~^ sr~^ for exam- 
ple) is also in use and the astro physical factor is ca n be found under 
this form in some articles (e.g., lEvans et al.ll2004l) . 



APPENDIX B: TOY MODEL FOR J (IN DSPHS) 

The volume of the dSph is not always fully encompassed in the 
integration solid angle, as sketched in Fig. IB II (vertical hatched 
region) so that a numerical integration is required in general. 
However, a reasonable approximation for estimating the depen- 
dence of J on the parameters of the problem, i.e. the distance to 



^' Some authors do ho wever expli c itly ex press the solid angle depen- 
dence in their units, e.g. jPieri et al who express J (<3? 

Cosmo in 

their notation) in GeV^cm~®kpc sr. This is completely equivalent to our 
JVfSjkpc^''' but for the unit numerical conversion factor 
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Figure Bl. Sketch of the integration regions contributing to the J factor: 
shown are the full integration region (vertical hatched) or a sub-region 
(cross-hatched) used for the toy calculations. The letter O shows the ob- 
server position, Oint is the integration angle, d is the distance of the dSph 
and -Rvir its virial radius. 



the dSph d, the integration angle aint, and the profile parameters 
Ps, Ts and 7), is to consider only the volume within the radius 



d X sin(aint) ~ d x aint, 



(Bl) 



where the approximation is valid for typical integration angles 
Qint 0.1°. This volume corresponds to the spherical cross- 
hatched region in Fig. IB II 

The toy model proposed below to calculate J allows us to 
cross-check the results of the numerical integration for both the 
smooth and sub-clump contribution. We find that the model is ac- 
curate enough up to a factor of 2 for 7 = and 7 > 0.5, so can be 
used for gross estimates of any signal from a DM clump. 



Bl For the smooth distribution 

About 90% of the clump luminosity is usually contained in a few 
rs, whatever the profile. The consequences are twofold. First, as 
can be read off Table [2] Vs/d <^ 1, so that the J factor amounts to a 
point like contribution 



J, 



point — like — 



47r /■min(rQ.^^,rs) 



d2 



r p {r)dr. 



(B2) 



Secondly, it means that Eq. ^ for the profile can be simplified into 
the approximate expression 



Papprox (r) 



Psat 
Ps X 





if r ^ rsat ; 

if r-sat < r- < rs 

otherwise. 



(B3) 



However, for all applications of our toy model, we will keep 7 < 
3/2, so that the saturation density above is never reached in the 
dSphs considered below. 

Various regimes The approximate formulae for J is obtained by 
combining Eqs. iB2i and ( |B3l l: 



47r /''"'"('■°i„t'''=) 



d2 



r Papprox(''-)dr. (B4) 



Using Eq. l IBlb . this leads to 

d2 '3-27 



J. 



approx — 79 



[min(rQi„j,rs 



13-27 



(B5) 



This formula gives satisfactory results for cuspy profiles (see be- 
low), but has to be modified in the following cases: 



• If rai,,! > rs, the integration region encompasses rs- The 
(1,3,7) profiles decrease faster than for r ~ rs hence in- 
tegrating the toy model up to rs is bound to overshoot the true 
result. We thus stop the integration at the radius r^ such that 

Ptruc(ra;) — Papprox (r^; ) /x, i.e. 

r, = rs • [x^^^'--'^ - 1] . 

Taking x — 2 gives a satisfactory fit to the full numerical calcula- 
tion (see below). 

• If '"Qint ~ '"s and 7 = 0, the integration can be performed 
analytically up to i?vir and is used instead. 

• If r^int and 7 = 0, the profile is constant, and integrating 
on the cross-hatched region (instead of the vertical hatched one, see 
Fig. IBU undershoots the true result. A better approximation is to 
integrate on a conic section. For the same reason as given for the 
first item, we replace rs by rx (with a; = 2) in the calculation of 
the cone volume. 

Resulting formula To summarise, the final toy-model formula 
proposed for the smooth contribution of the dSph is: 



Hoy 



where 



47rp| 
d2 



X < 



„27 



i(r 



X3-27 



3-27 

ri_ • rs 



if 7 > 0; 

if 7 =0, rc,i„t >r^; 

if7 =0, r<,i„^ <rx- 
(B6) 



= aint • d, 



1/(3-7) 



1], 



I{x) 



-rt{rt + 5rsX + 102;')/(30(rs -f xY). 



(B7) 



Toy model vs numerical integration Finally, we check the va- 
lidity of this toy model by confronting it with the full numerical 
integration. Various inner slope 7 of the profile are considered as 
provided in Table [T] Defining the critical distance dcrit for which 
the dSph is fully encompassed by the integration region, i.e., 

, _ rs 

Q-crit — 

aint 

we find dcrit ~ 50 kpc and 500 kpc for rs = 0.1 and 1 kpc respec- 
tively (the integration range is Ofint = 0.1°). If r^ is used instead of 
rs, this distance is even smaller. This allows us to test the toy model 
for the two regimes. The result is shown in Fig. IB2I The symbols 
show the full numerical integration while the lines show the toy- 
model calculations. For profiles steeper than 0.5, the agreement is 
better than a factor of 2 for all distances. For flatter profiles, the toy 
model only gives results within an order of magnitude. However, 
for 7 = 0, the fix applied to the toy-model allows to regain the 
correct results within a factor of 2. 

Hence, given the current uncertainties on the profiles, the set 
of formulae jBll and l |B7t can safely be used for quick inspection 
of the J value of any profile with an inner slope 7 of 0, or greater 
than 0.5. 



B2 For the sub-clumps 

The influence of DM sub-structures on the 7-ray production has 
been widely discussed in the literature. These sub-structures may 
enhance the detectability by boosting the 7-ray signal. In this ap- 
pendix, we give an analytical estimation of the effect of sub-clumps 
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Figure B2. Toy-model calculation (lines) vs full numerical integration 
(symbols) of J as a function of the distance to the dSph. The integration 
angle is fixed to ajnt = 0.1° and the (1, 3, 7) profiles are taken to vary 
from 7 = to 7 = 1.45. For each model, pa is calculated such as to 
provide M300 = lO'^ Mq. Top: dSphs for which rs = 0.1 kpc. Bottom: 
dSphs for which r's = 1 kpc. . 



in dSph spheroidal galaxies, in the same spirit as the toy model de- 
veloped in the previous section for the smooth component. For sim- 
plicity, we restrict ourselves to one cored (a, /3, 7) — (1, 3, 0) and 
one cusped (1, 3, 1) profile. T o characterise t he clu mp distribution, 
we use the formalism given in lLavalle et al ] (l2008h . 

Sub-structure distribution The clump spatial distribution is as- 
sumed to follow the dSph DM profile, namely 



dP(r) I r 

^ oc 

dV 



1 + 



(B8) 



The mass distribution of the clumps is taken to be independent of 
the spatial distribution and takes the usual form, 

dP 



dM 



= AM-" 



(B9) 



with Af £ [Afniin, Mn-iax] and a ~ 1.9 from cosmological N-body 
simulations (A is the normalisation constant for dP/d.M to be a 
probability). 

Clump luminosity Defining Li the intrinsic luminosity of the sub- 
clump i to be 



p'dV, 



(BIO) 



the astrophysical contribution to the 7-ray flux from the sub- 
structures of the dSph is 



J clun 



d2 



(BU) 



where N'^ is the number of clumps contained within the integration 
angle a and d is the distance of dSph. The luminosity depends only 
on the mass of the clump, once a concentration-mas s (cvir — Afvir) 
relationship is chosen (see, e.g. jLavalle et al .I2OO8I and references 
therein), so that Li = L{Mi). Moving to the continuous limit, 
Eq. jBllb reads 



r _ _j_ ivrcl 

f-'clumps JO 



LiM)§dM 



(B12) 



Fitting the results from lLavalle et al. I( l2008h . the intrinsic luminos- 



it>|_| varies almost linearly with the mass of the clump, as 
L'^^'*(M) = 1.17 X 10* {M/Mq)°-^^ Ml kpc~^ 



so we have 



J clun 



N^^A /1.17 X 10* 



d2 



1.91 - a 



M, 



M 



(B13) 



(B14) 



Number of clumps The fraction F of clumps in the spherical in- 
tegration region rcj^t « Ointd (cross-hatched region in Fig. ID IS 
given by 



F = 



dV 



(B15) 



where A'^^'ot is the total number of clumps within the dSph. Upon 



integration and defining x-mt = ra 
becomes: 



j/rs and x^ir = Rvii/rs this 



Fc. 



and 



4:Xc+S 3 

J + ln(a:c, + 1) - - 



2(a;c + 1) 
4a;vir + 3 



2(,,, + l)2+l^(-vir + l)-- 



-fcusp — 



1 



(Xa + 1, 
1 



1 



+ ln(a:„ + 1) - 1 



+ ln(a;vir + 1) - 1 



(B16) 



for (1,3,0), 



(B17) 



for NFW. 



Some care is necessary when evaluating the number of clumps 
N'^^ — Fx Ntot in the integration region. Whatever the profile, 
most of the clumps are located within rs so when ra^^^ > r^, the 
spherical integration region of our toy model (cross-hatched region 
in Fig. IB lb is a good enough approximation, and Eq. l lB16t and 
l |B17| ) hold. However, if ra-^^^ < Vs then the remainder of the in- 
tersecting cone (vertically hatched region in Fig. IB It could amount 
to a significant contribution to the number of clumps. Cuspy dis- 
tributions should only be marginally affected given their high cen- 
tral concentration. However, this effect may be important for cored 
profiles. Whenever ra^^^ < r^, as for the smooth contribution, 
Eq. l |B16l ) is therefore multiplied by the ratio of the intersecting 

In this toy model, we limit ourselves to the NFW profil es for the sub- 
clump s in the dSph, and a Cvir — Afvir relation taken from lBullock et all 
i200ll) . 
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Figure B3. Toy-model calculation (lines) vs full numerical integration 
(symbols) of J as a function of the distance to the dSph. The integration 
angle is fixed to a-^^^ = 0.1° and the two (1, 3, 7) sub-clump spatial distri- 
bution are 7 = a nd 7 = 1 (their inner profile is a NFW with a Cvir — A/vir 
relation taken from lBuUock et al1j20Qlb ). The calculations assume the frac- 
tion of DM in sub-clumps to be / = 50% of the total mass of the dSphs, 
where the smooth profile is taken as in Fig. HI] Top: Ta = 0.1 kpc. Bot- 
tom: rg = 1 kpc. 



cone volume to the integration sphere volume, in order to account 
for that effect. 

If the mass of the dSph is A/vir and assuming a fraction / of 
this mass is in the form of clumps, one gets using Eq. ( |B9l l 

Nti = / ^A/vi. {Mf-^: - Mi^:)-' . 

Resulting formulae Adding all ingredients together, the contribu- 
tion of the sub-structures to the flux is 



Jclumps = 1.17 X 10 



8 -^core/cusp f d 



M, 



1.91-a _ 71 *-1.91-a 
max ^^^min 



f1-~a Ttjl — a 
max ~~ min 



1.91 - a 



(B18) 



Toy model vs numerical integration The comparison between 
the two is shown in Fig. IB3I The symbols show the full numeri- 
cal integration while the lines show the toy-model calculations. For 
ra = 100 pc, the agreement is better than a factor of 2 for all dis- 
tances. For rg = 1 kpc, the toy model only gives results correct to 
within a factor of 4 for 7 = 1. 

Hence, given the current uncertainties on the profiles. 
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Figure CI. Jsm(S = 0) as a function of the distance to the dSph for three 
profiles 7 and three values of . The corresponding values for ps are given 
in Table[T] 



Eq. l lB19t can be used for quick inspection of the J value for the 
sub-clump contribution. 



APPENDIX C: DISTANCE AND INTEGRATION ANGLE 
DEPENDENCE ON J FOR GENERIC DSPHS 

This Appendix completes the study of the J-factor dependences 
started in Section l23l All the plots and discussions below rely on 
the generic profiles given in Table[T] and the sub-structure reference 
configuration given in Section l2.2.3l 



CI Distance dependence J(d) 

Fig. lCll shows Jsm as a function of the distance to the dSph (we as- 
sume Qint = 0.1° here and that we are pointing towards the dSph 
centre, i.e. Q = 0). As we have checked earlier, the sub-clump con- 
tribution for the reference model at 6 = is always sub-dominant, 
so for clarity only Jsm is displayed (/ = 0) in the figure. 

If the angular size of the signal is smaller than the integration 
angle, the distance dependence is expected to be Jsm oc d This 
is the case for 7 = 1.5 for any value of Ts (hollow squares curves). 
Actually, the three curves follow the point-like source toy formula 
( |B6t appropriate for steep 7, i.e. 



j{e = 0) oc Ps X 



d? 



(Cl) 



However, when the angular size of the emitting region becomes 
larger than the integration angle, the above relationship fails. As 
most of the flux is emitted within js, this happens for a critical 
distance 



Cfint 



(C2) 



For rs = 0.1 kpc, this comsponds to dcrit ~ 60 kpc (see the full 
circles dashed curve for 7 = 0). Having a dSph closer than this crit- 
ical distance does not increase further the signal (see, e.g., the solid 
and dotted full circles curves for 7 = and > 0.5 kpc). In the 
latter case, taking a larger integration region is not always the best 
strategy as, from an experimental point of view, a larger integration 
region increases not only the signal but also the background. In this 
case, the gain in sensitivity from having a dSph close by is not as 
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Figure C2. J x (d/100 kpc)^ as a function of Oint x (d/100 kpc) for a 
generic dSph with rs=0.5 kpc: smooth (thick dashed lines) and sub-clumps 
(thin dotted lines). With this rescaling, the case d = 10 kpc (stars) super- 
imposes on the case d = 100 kpc (empty and full circles). 



important as what might naively be expected from the point-like 
approximation (see Section[3). 



C2 Integration angle dependence J(Qint) 

We recall that /^q dJl = J^^ dfiint Jq"'"' sin(aint)daint, where 
AQ = 27r(l — cos(Qint)), so that the J-factor from Eq. (|5} can be 
rewritten in the symbolic notation 



(C3) 



with 



and 



[/3i„tl 



(C4) 



Jo 

Jo 

-f'lAnt.aint] = sin(ai„t) / J'[r(/, /3int, aint)] d/. (C5) 
Jo 

For small integration angles and the case of a flat enough profile, 
the integrand in Eqs. iC4\ and dCSt does not vary much with aint, 
so that for the smooth ( J^ = P^) and the mean sub-clumps = p), 
we have 



J SI 



and (Jsubci) oc af„ 



(C6) 



Fig. IC2I shows the integration angle dependence for the 
smooth (1 — /)^Jsm (dashed lines) and the sub-clump mean 
(Jsubci) (dotted lines) contributions. (The pointing direction is to- 
wards the dSph centre.) For 7 = (solid black circles), the af^t 
scaling holds up to af^f ~ 3° if d = 10 kpc (as given by Eq. lC2t . 
A plateau is reached when the entire emitting region of the dSph is 
encompassed (i.e. for a few Vs/d). For 7 = 1 (blue empty circles), 
the curves are slightly more difficult to interpret, as the profile is 
not steep enough for it to be considered fully point-like (and thus 
'independent' of Qint) given the integration angles considered. 



Table Dl. Maximum boost and transition regime, i.e. (ctintd-) B—l in 
deg kpc, for which B = 1, for various smooth/sub-clump parameters for 
three inner slope 7 (for the smooth). 



Config.t 



7 = 



7 = 0.5 



Bn 



reference-^ 


1.9 1 19 


2.2 


21 


2.0 1 30 


[global parameters] 












Q = 1 


1.0 


40 


1.3 


60 


1.6 1 160 


13 = 5 


2.3 


11 


2.0 


18 


1.3 1 36 


Rvlr = 6 kpc 


3.0 


15 


3.5 


20 


2.9 29 


M300 = 2 ■ 10'' Mq 


1.3 


66 


1.4 


52 


1.3 64 


[sub-clump parameters] 












dP/dV =Einasto* 


1.4 




1.7 




1.7 1 22 


a = 1.7 


1.3 


62 


1.5 


50 


1.3 61 


a = 2.0 


2.8 


0.2 


3.4 


8 


2.9 16 


M^in = IMq 


1.5 


43 


1.7 


37 


1.5 47 


Afmax = W^Mq 


2.4 


4 


2.8 


14 


2.5 22 


/ = 0.5 


3.4 


10 


4.2 


16 


3.5 25 


Psubcl =Einasto 


8.7 


0.05 


10.6 


0.35 


9.0 4 


Cvii- X 2 


7.6 


0.06 


9.3 1 0.4 


7.9 1 4.5 



t All parameters are as for reference, except those quoted, 
t Reference configuration (M300 = 10^ Mq): 

■ pBui with (a, /3, 7) = (1, 3,7) and dP/dV oc Psm \ 

■ i?vir = 3 kpc and rs = 1 kpc (for pam and dP/dV); 

■ dP/dM = Af-° (a = 1.9), and M^uh £ [10" 

■ / = 0.2, ps^^bcl =NFW, and c^j, - M vi,=B01 
* Einasto parameters taken from lMerritt et al.l l l2006l) . 



Finally, the rescaling used in Fig. lC2l implies: 

Jdl(aint) = Jd2 (aint|) X (I 



APPENDIX D: COMPLEMENTARY STUDY OF THE 
BOOST FACTOR 



(C7) 



In Section [2.3.21 we concluded that the boost could not be larger 
than a factor of 2 for all configurations where the sub-clump spa- 
tial distribution follows that of the smooth halo in the dSph. The 
calculations were also made for a 'reference' configuration of the 
sub-clumps. However, the boost can be smaller (or larger) when the 
latter parameters are varied. 

In Tab. ID II we systematically vary all the parameters entering 
the calculation in order to compare with the reference model case. 
The two quantities of importance are the maximum boost possible 
(which is obtained when aint fully encompasses the clump), and 
the transition point aintd for which the boost equals 1 (the min- 
imum value is always given by (1 — f )^)- The reference results 
correspond to the numbers obtained from the dotted lines in Fig. [5] 
i.e. for Ts = 1 kpc. Note that most of the values for Bmax in the 
Table would be close to unity if rs = 0.1 kpc were to be selected. 



The dependence can be understood by means of the toy model formulae 
)B6> and )B7t . For Oint < "^hii*. we have 

'^[7>0.51 <^ X (aintrf)^"^^ . 

For 7 = 1 (empty blue circles), J is then expected to scale linearly with 



Qint, which is observed for the smooth (dashed blue line), and to some 
extent for the sub-clump contribution (dotted blue line). However, for the 
latter, the transition region (around rs) falls from a slope a = 1 towards an 
outer slope /3 = 3 (instead of falling from = 1 to = 6. Hence, for 
Qint > crfnt' ' ''^^ sub-clump contribution continues to build up gradually. 
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Dl Varying the [global parameters] 

The four lines under '[global parameters]' keeps the recipe of 
dP/dV psm^ but some previously fixed parameters are now var- 
ied. The trend is that a sharper transition zone (larger a), a larger 
radius of the dSph, or a smaller mass imply a larger Bmax- The im- 
pact of the outer slope /3 depends on the value of the inner slope 7. 
However, the maximum boost factor reached for these parameters 
is never larger than ~ 3. The typical transition value lies around 
20° kpc, which corresponds, for a dSphs located 100 kpc away, 
to an integration angle of 0.2°. Hence, for all these configuration, 
large integration angle should be preferred (this is even worse for 
closer dSph). 

D2 Varying the [sub-clump parameters] 

The remaining lines under [sub-clump parameters] show the im- 
pact of the choice of the distribution of sub-clumps, the mass dis- 
tribution parameters (minimal mass and maximal mass of the sub- 
clumps, slope a of dP/dM), and the density profile of the sub- 
clumps. Relaxmgjte condition dP/dV oc psm has no major im- 
pact. In lSpringel eTal] ilOOA a simple Einasto profile with univer- 
sal parameters was found to fit all halos (from the Aquarius sim- 
ulation) independently of the halo mass. For that spec i fic cas e, we 
use the values found for the Galaxy in iMerritt et all ( l2006h . The 
Einasto profile is steeper than 7 = but it decreases logarithmi- 
cally inwards. Only for 7 > 1 (for the smooth component) such 
a model is able to marginally increase the maximum boost w.r.t. 
the reference model (instead of decreasing it), which is not unex- 
pectedS Varying the mass distribution slope a is understood as 
follows: for a ~ 1.9, all decades in mass contribute about the same 
amount. When a is decreased, the less massive sub-halos dominate, 
whereas for a > 1.9, t he most massive su b-halos dominate the lu- 
minosity (e.g. Fig. 4 o f iLavalle et al.l2008l) . This has to be balanced 
by the fact that the fraction of DM going into sub-clumps remain 
the same (/ — 0.2), regardless of the value of a, so that the total 
number of clumps in a mass decade also changes. The net result 
is a smaller boost when a is decreased, and a larger boost from 
the more massive sub-structure when a is increased. In a similar 
way (a is now fixed to 1.9 again), the mass also impact on B, but 
in a marginal way. The only sizeable impact comes from varying 
the fraction of mass into clumps, the sub-clump profile or the con- 
centration of sub-clumps. In the first case, when / increases, the 
smooth signal decreases by (1 — f)^ whereas the sub-clump signal 
increases as /. Even if / is incre ased up to 50%, whic h is very un- 
likely (recent simulations such as lSpringel et al.l ( f2008h tend to give 
an upper limit of / < 10%) this gives only a mild enhancement. In 
the second configuration, the NEW profile for the sub-clumps is re- 
placed by an Einasto one. Despite its logarithmic slope decreasing 
faster than the NEW slope 7 = 1 below some critical radius, the 
latter profile is known to give slightly more signal than the NEW 
one (pEinasto(r) > Pnfw('') for a region that matters for the J 
calculation). This results in a boost close to 10, regardless of the 
dSph's smooth profile. Finally, we recall that the BOl Cvir — M^h 
relation is used to calculate the value of the scale parameter for 
any sub-clump mass. In the last configuration, the concentration 

For smaller 7, the smooth distribution, in that case, is flatter than the 
sub-clump one, so that the boost is larger than one for small ajnt and the 
transition where B = 1 is ill-defined. However, such a configuration is 
highly unlikely as it is exactly the opposite of what is observed in all N- 
body simulations. 
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Figure El. 80% containment radius (9go) of PSF-convolved DM annihila- 
tion halo models versus ago ■ 

parameter is simply multiplied by a factor of 2, which is probably 
not realistic. Again, the same boost of ~ 10 is observed. Accord- 
ingly, for these last two cases, the transition angle is reduced, and 
corresponds to Qint < 0.01° (for a dSphs at 100 kpc). 

To conclude, although boosts by as large as a factor of 10 can 
be obtained through suitable combinations of parameters, most of 
these combinations are unlikely and require the signal to be inte- 
grated on large angles. 



APPENDIX E: IMPACT OF THE PSF OF THE 
INSTRUMENT 

Fig. lEll shows the impact of the instrument angular resolution on 
the 80% containement radius for J (for the generic dSphs studied 
in Sec. [3). The solid line corresponds to the quadrature approxi- 
mation given by Eq. J15t . whereas the symbols correspond to the 
convolved PSF*halo profile. The PSF is described by the sum of 
two Gaussians and is a scaled (factor two improved) version of the 
PSF appropriate for H.E.S.S. at 200 GeV. Calculated halos for a 
range of a, /?, 7 models consistent with the stellar kinematics of 
the classical dSphs are shown as gray squares. The quadrature sum 
approximation used in this work is shown as a solid line. 



APPENDIX F: CONFIDENCE LEVELS AND PRIORS 

In this Appendix, we describe how confidence intervals for the 
quantities such as p(r) or J are chosen. 

Fl Sensitivity of the result to the choice of prior 

In the Bayesian approach, the PDF of a parameter x is given by 
the product of the MCMC output PDF P{x) and the prior p{x). 
The resulting PDF is therefore subjective, since it depends on the 
adoption of a prior. However, whenever the latter are not strongly 
dependent on x, or if ^(2;) falls in a range where p{x) does not 
strongly varies, the PDF of the parameter becomes insensitive to 
the prior. This happens for instance if the data give tight constraints 
on the parameters. 

In our MCMC analysis, we assumed a flat prior for all our 
halo parameters, as there is no observationally motivated reason 
for doing otherwise. Note, however, that flat priors on the model 
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parameters do not necessarily translate into flat priors on quantities 
derived from those parameters. Specifically, the flat priors on our 
model parameters imply a non-flat prior on the DM density (and 
also on its logarithm) at a given radius, and hence a non-flat prior 
on J. In principle it is possible to choose a combination of priors 
for the parameters that would translate into flat priors on p{r), but 
we have not done so here. The general impact of such choices, and 
the methodology to study the prior-dependent res ults, has been dis- 
cussed in the context of cosmological studies bv lValkenburg et al] 
( I2OO8I) . In this study, we only use a flat prior on the parameters (or 
on the log for Vs and ps). The test with artificial data demonstrate 
that our reconstructed p and J values are sound. 

F2 Confidence intervals for p(r) and cross-cliecks 

F2.1 Definition 

Confidence intervals (CI), associated with a confidence level 
x% (CL), are constructed from the PDF. The asymmetric interval 
Ax = [Ox , 9x] such as 

CL(a;) = [ V[6)Ae= 1 - 7, 

J A,. 

defines the 1 — 7 confidence level (CL), along with the CI of the 
parameter 6. We rely on two standard practises for the CI selection. 
The first one (used only in this Appendix) is to fix 6~ to be the 
lowest value of the PDF. The CLs correspond then to quantiles. 
This is useful for CI selection of values, to ensure that the best- 
fit value of a model (i.e. t he lowest x^) falls in the CI (see, e.g., 
Fig. 7 of IPutze et al.ll2009h . In the second approach, the CI, i.e 6^ 
(resp. 9x), is found by starting from the median g™^'* of the PDF 
and decrease (resp. increase) 6x until we get x-%/2 of the integral 
of the PDF. This approach ensures that the median value of the 
parameters falls in the CI, any asymmetry in the CI illustrating the 
departure from a Gaussian PDF: this is the one used thoughout the 
paper. 
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Figure Fl. Projected distribution of logj^Q(p) along with the value of sev- 
eral other estimators for the MCMC analysis of Draco. In this box projec- 
tion, the larger the box, the most likely the probability of log]^Q(p). For 
instance, on the top panel, for logig(r) = —1.5 the probability density 
function of logj^Q(p) is distributed in the range [8 — 10] and peaks around 
9.5). 

The CLs obtained from it are larger than the previous one. In the 
above paper, the uncertainties were small even with that method, 
so that was not an issue. However, in this study, this makes a huge 
difference in the resulting value CL of J. 

In order to check which approach was the correct one, we 
bootstrapped the Draco kinematic data and calculate from the col- 
lection of p(r) from each bootstrap sample the median value and 
the uncertainty. The first approach, where the CLs are directly cal- 
culated from the full set of MCMC samples was in agreement with 
the bootstrap approach, meaning that the second one biases the re- 
sults toward too large uncertainties. The results of the paper rely 
thus on the first and correct approach. 




F2.2 Comparison of several choices for the PDF of p(r) 

Fig. lFll shows the projection for each r of the PDF calculated from 
the output MCMC file. To do so, p(r) is calculated for each en- 
try of the thinned chains and then stored as an histogram. This re- 
sults in 'boxes': the larger the box, the more likely the value of 
p(r). From this distribution, we can calculate the median (thick 
solid black line), the most probable value (thick dotted black line). 
The thick solid red line correspond to the model having the small- 
est value. We see that the latter differ from the median one for 
this dSph, though they can be close for other dSph in our sample. 
In this paper, as our analysis is based on the Bayesian approach, we 
disregard the best-fit model and only retain the median value. 

In the first approach, the 68% and 95% CLs are calculated 
from the distribution pr (at each r) They are shown as dashed and 
dotted thick black lines. Note that none of all the above lines corre- 
sponds to a physical configuration of p{r). 

A second approach is to construct the 68% CLs from a sam- 
pling of the (still) correlated parameters. This is achieved by using 
all sets of parameters {0}x%Ch ~ {^j}i=i ■ p, for which x^{Oi) 
falls in the 68% CL of the x^ PDF (see above). Once these sets are 
found, we calculate p(r) for each of them, and keep the maximum 
and minimum values for each position r. This defines envelopes 
of p(r) (CIs are found for each r). This is s hown as dot t ed an d 
dashed red lines. Such an approach was used in lPutze et al. 



APPENDIX G: ARTIFICIAL DATA SETS: VALIDATION 
OF THE MCMC ANALYSIS 

In this section, we examine the reliability of the Jeans/MCMC anal- 
ysis by applying it to artificial data sets of 1000 stellar positions 
and velocities drawn directly from distribution functions with con- 
stant velocity anisotropy. We assume the form L~^'''""=°/(e) for 
the distribution functions, where the (constant) velocity anisotropy 
is given by /?aniso = 1 — cr'^/a'^, with and a'^rnr being the 
second moments of the velocity distribution in the radial and tan- 
gential directions, respectively. The function /(e) is an unspeci- 
fied function of energy e which we determine numerically using an 
Abel inversion once t he halo model and stellar dens ity are speci- 
fied (ICuddefordlll99lh . We used the same models in I Walker et al.l 
1 I20T l|), but we present here a more general study. The set of artifi- 
cial data covers a grid of models with 7 = 0.1, 0.5, 1.0, rh/rs = 
0.1, 0.5, 1.0 and /3 = 3.1. For each halo model, we assume /3aniso 
values of (isotropic), 0.25 (radial) or —0.75 (tangential): the 
/3aniso values for the anisotropic models are chosen to give mod- 
els with roughly equivalent levels of anisotropy (in terms of the 
ratios of the velocity dispersions in the radial and tangential di- 
rections). We also generate a grid of models with a steeper in- 
ner slope 7 = 1.5 and /3 = 4.0. In all cases, the haloes contain 
~ IO^Mq within 300pc. We mimic the effects of observational er- 
rors by adding Gaussian noise with a dispersion of 2 km s^ ^ to each 
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individual stellar velocity generated from the distribution function. 
The reconstruction depends on the choice of the prior 7prior, and 
this effect is explored in the two sections below. 



Gl Prior: ^ 7prior ^ 1 versus < 7prior ^ 2 

We start with the free 7prior analysis (see Sec. 15. lb based on two 
different priors. Top panels of Fig. iGTI show the ratio of the recon- 
structed median profile to the true profile. There is no significant 
differences for p(r > 1 kpc) when using the prior ^ 7prior ^ 2 
(top right) or ^ 7prior ^ 1 (top left); at large radii, the pro- 
file does not depend any longer on the 7 parameter. However, it is 
striking to see that restricting the prior to ^ 7prior ^ 1 greatly 
improves the determination of the inner regions for the profile, re- 
gardless of the value of 7truo. Even for 7truc = 1-5 (green curves), 
using an incomct prior does not degrade to much the reconstruc- 
tion of the profile. 

This results is further emphasised when looking at J. The bot- 
tom panels of Fig. lGll are plotted with the same scale to emphasise 
the difference. As J integrates over the inner parts of the profile, 
the median MCMC value can strongly differ from the true value 
for cuspy profiles. This difference can reach up to 5 orders of mag- 
nitude (over the whole range of aint) for 7truo ^ 0.5 when using 
the prior 7prior ^ 2 . The prior 7prior ^ 1 does generally 
better, and accordingly, the confidence intervals are much smaller 
than for the other prior (for any integration angle). 

The behaviour of the 7truo = 1-5 case is unexpected. Using 
the prior ^ 7prior ^ 1 does better than the other one for any 
integration angle. Indeed, even if the reconstructed median value 
is shifted by a factor of 10, its CLs correctly encompass the true 
value. It does better than the ^ 7prior ^ 1 prior, which correctly 
provides CLs (that bracket the true value), but which are completely 
useless as these CLs can vary on ~ 8 orders of magnitude. 



G2 Strong prior: 7prior Axed 

In Fig. IG2I below, we use a prior 7prior = for models having 
7true = 0, a prior 7prior = 0.5 for models having 7truc ~ 0.5, etc. 

A comparison of Figs. IGll (using 7prior ^ 1 or ^ 
7prior ^ 2) and IG2I (fixed 7prior) shows that the latter prior only 
slightly improves the precision of the J-factor reconstruction for 
7true = 0, 7truo = 0.5, and 7truc = 1. However, if 7truc = 1-5 
(green curves), although the corresponding J-factor is now better 
reconstructed than when using the prior ^ 7prior ^ 2 (Figs lGll 
top panel), it is suiprisingly less reliable than the strongly biased 
^ ^ 1 prior. 

The main conclusion is that the knowledge of 7truG does not 
help providing tighter constraints on J: the uncertainty remains 
a factor of a few, except when the inner profile is really cuspy 
(7truc = 1.5), in which case it becomes strongly biased/unreliable. 



APPENDIX H: OTHER RECONSTRUCTION 'BIASES' ON 
THE J-FACTOR 

In this Appendix, the MCMC analysis is performed based on the 
prior ^ 7prior ^ 1, for which the analysis is found to be the most 
robust (see previous Appendix). 
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Figure G2. Fixed 7prior MCMC analysis. Top: p{r). Bottom: J(aint). 
HI Impact of tlie binning of the stars 

Figure IHll shows the impact of using different binnings in the 
MCMC analysis. The left panel shows the reconstructed (median) 
value of the velocity dispersion as a function of the logarithm of r 
(to emphasise the differences at small radii), for a binning used in 
this paper (black; where each of \fN bins has VN member stars, 
where A'' is the total number of members), a binning with two times 
(red) and four times (blue) fewer bins. For Fornax and Sculptor, the 
profiles are insensitive to the binning chosen, so that the reconstruc- 
tion of the J values median and 68% CLs (right panel) is robust. 
For other dSphs, either the adjusted velocity dispersion profile is 
affected at small radii, or at large radii. In the latter case, the J cal- 
culation should not be affected, as the outer part does not contribute 
much to the annihilation signal. In the former case, a deviation even 
at small radii can affect the associated J by a factor of a few. The 
exact impact depends on the integration angle, the distance to the 
dSph (which corresponds to a given radius), and the 'cuspiness' of 
the reconstructed profile (the J value of a core profile will be less 
sensitive to differences in the inner parts than would be a cuspy pro- 
file). For instance, Draco and Leol both have a 2 km/s uncertainty 
below 100 pc, but Draco is three times closer than Leol: their J for 
a given aint have different behaviours (right panel). The strongest 
impact is for Leol that have the fewest data. The flatness of the J 
curve seems to indicate a cuspy profile (all the signal in the very in- 
ner parts), which we know are the least well reconstructed ones (see 
Appendix IGU . Leol is thus the most sensitive dSph to the binning, 
for which a balance between a sufficient coverage over r and small 
error bars cannot be achieved. The ultra-faints dSphs are expected 



Dark matter annihilation in dwarf spheroidal galaxies and ■~^-ray observatories: I. Classical dSphs 29 




Figure Gl. Ratio of the MCMC profile to the true profile. The lines are colour-coded with respect to the value of the true inner slope 7truc of the artificial 
data. Top panels: ratio of the median p(r). The two vertical gray dashed lines correspond to the typical range within which the artificial data bin are taken. 
Bottom panels: ratio of J(oint) for the artificial dSphs located at 100 kpc. Left panels: MCMC analysis with the prior 7piior ^ 2. Right panels: the 
prior is ^ 7prior ^ 1- 



to have even fewer stars, so that their J calculation is expected to 
be even more uncertain. 

Overall, the choice of the binning can produce an additional 
bias of a few on the J reconstruction. This is an extra uncertainty 
factor that makes Fornax and Sculptor the more robust targets with 
respect to their annihilation signal. Surveys in the inner parts and 
outer parts of Carina, Draco, Sextans, Leo I, Leo II and Ursa Minor 
are desired to get rid of this binning bias. 



H2 Impact of the choice of the light profile 

Figure llI2l shows the various median values and 68% CIs of J when 
changing the assumptions made on the light profile. The black lines 
labeled 'physical' correspond to the Plummer model used for the 
main analysis (see Eg. I21t: the red lines labeled 'unphysical' are 
also Plummer, but the physical constraint given by Eq. ( 125) is re- 
laxed; the blue lines and green lines correspond respectively to a 
light profile modeled with an (a, /3, 7) profile in order to get a 
steeper outer slope (2, 6, 0) or a steeper inner slope (2, 5, 1) with 
respect to the Plummer profile. Regardless of the light profile used, 
we recover similar critical angles for which J is the most con- 
strained. The impact on the J value is strongest for the least-well 
measured profiles (Leo 2 and Sextans), but is contained within the 
95% CI and marginally within the 68% CL. 
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data). Right: coiTesponding median values and 68% CI for J(aint). 
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